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rilHE following lectures were collected from manuscripts left 
-*- by the late Professor J. C. Adams, and are now reprinted 
without change from his Collected Scientific Papers, Vol. ii., 
pp. l-r-84. 

It was thought that the wide interest attaching to the 
lunar problem reached many besides the professed astronomer, 
and would justify a separate publication of this short work. 

It is known that Adams contemplated the publication of 
some such essay himself, and it must be a matter of regret to 
all that he never did so. No pains have been spared to 
present the material properly, but it is unavoidable that it 
should apgear from the hands of an editor in a less perfect 
form than if the author had issued it himself. 

Yet, allowing for this disadvantage, I think those best 
qualified to judge will consider this work fally worthy of 
Adams's great name. Of current elementary theories it may 
be said that they leave oflf where the difficulties of the subject 
begin, that is to say, where the various cases of slow con- 
vergence have been exposed, but not dealt with. It is perhaps 
not too much to say that these lectures carry us to the point 
where such difficulties end, in an adequate evaluation of all 
the chief constants. They leave the problem effectively solved 
and not merely stated, and shew the path clear for the for- 
mation of a detailed theory, if that is desired. 

R. A. SAMPSON. 

Durham. 

8 October, 1900. 
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2 LECTURES ON THE LUNAR THEORY. 

a more adequate knowledge than we now possess of the relative 
motion of Three Bodies. So far, with the notable exception 
of Dr G. W. Hill and those that have followed him, we have 
seen comparatively little effort in this direction. 

This was the cardinal feature of Adams's plan, and his 
lectures shew the methods he had gradually elaborated to 
accomplish it. They separate the inequalities from one another 
as far as possible, and are content with indicating the manner 
in which these separate inequalities afterwards combine. To 
shew that, with so slight an apparatus and within so small 
a compass, the result is no mere sketch, we need but set 
side by side the coefficients of longitude found in these 
Lectures and the corresponding terms in Delaunay's Thiorie, 

Adams. Delaunay. 
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For those to whom the difficulties of the Lunar Theory 
are known, these numbers need no comment. 

No Manuscript exists of Lecture I. It is taken substantially 
from my own notes of 1889.] 

♦ With Delaunay'8 value of the Sun's Parallax, viz. 8"-76. 



LECTURE I. 

HISTORICAL SKETCH. 

[The Lunar Theory may be said to have had its commence- 
ment with Newton. Many irregularities in the Moon's motion 
were known before his time, but it was he that first explained 
the cause of those irregularities and calculated their amounts 
from theoiy. 

Of the inequalities which are due to the action of the 
Sun, the first, — which is called the Evection, — was discovered 
by Ptolemy, who lived at Alexandria in the first half of the 
second century of our era, under the reigns of Hadrian and 
Antoninus Pius. At a very early period the relative distance 
of the Moon at different times could be told from the angle 
it subtended, and its orbit could thus be mapped out. By 
such means Ptolemy found that its form was not the same 
from month to month, and that the longer axis moved con- 
tinually though not uniformly in one direction. He represented 
this change by a motion of the centre of the ellipse, as we 
would put it, in an epicycle round the focus, obtaining thus 
a variable motion for the longer axis and a variable eccentricity. 

The representation of position by means of epicycles is 
intimately related to the modem method of developing the 
coordinates in harmonic series; thus if we have 

a? = ill cos (^1^ + OLj) -f J.2COS {n^t 4- Oa) + ... 

y = Ai sin (riit + oli) + A2 sin (nj; -h 02) + . . • 

the motion of the point (a?, y) is that on a circle of radius ili 
with angular velocity n^, around. a centre which moves on a 

1—2 



4 LECTURES ON THE LUNAR THEORY. [LECT. 

circle of radius A^ with angular velocity n^, and so on; and 
if, more generally, we have 

a? = ill cos (M + tti) + ... 

y = Bi sin (nj; 4- aO + ... 

we may reduce this case to the former by rewriting 

a? = 2 (-^1 + Bi) cos (riit + tti) 4- 5 (Ai - B^) cos (- Wi^ - aO + ... , 

y = 2 (A^ + Bj) 8m(nit + aO + 2 (-^1 - ^1) sin (- n^^ - ai) + ... . 

Probably we have here the reason why circular motions 
and epicycles were first employed. 

Tycho Brahe (1546 — 1601) discovered the existence of 
another inequality in the Moon's Longitude quite different 
from the Elliptic Inequality and the Evection. He found it 
bore reference to the position of the Sun with regard to the 
Moon ; so that when the Sun and the Moon were in conjunction 
or opposition or quadratures the position of the Moon was 
quite well represented by the existing theory, but from con- 
junction to the quadrature following, her position was more 
advanced than the place assigned to it, reaching a maximum 
of some 35' about half-way; and in the second quadrant it 
was just as much behind. This inequality he called the 
Variation; it was the first that Newton accounted for theo- 
retically, and if we were to suppose the Moon and Sun to 
move, except for mutual disturbance, in pure circles in the 
same plane, it is the only one that would present itself. 

The next significant step was made by Hon-ox (1619 — 1641) 
who represented the Evection geometrically by motion in a 
variable ellipse, and gave very approximately the law of varia- 
tion of the eccentricity and the motion of the apse. He 
supposed the focus of the orbit to move in an epicycle about 
its mean place. 

Newton's Principia did not profess to be and was not 
intended for a complete exposition of the Lunar Theory. It 
was fi^gmentary; its object was to shew that the more 
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prominent in'egularities admitted of explanation on his newly 
discovered theory of universal gravitation. He explained the 
Variation completely, and traced its eflfects in Radius Vector 
as well as in Longitude; and he also saw clearly that the 
^change of eccentricity and motion of the apse that constitute 
the Evection could be explained on his principles, but he did 
not give the investigation in the Prindpia, even to the extent 
to which he had actually carried it. The approximations are 
more difficult in this case than in that of the Variation, and 
require to be carried further in order to furnish results of 
the same accuracy as had already been obtained by Horrox 
from observation. He was more successful in deaKng with 
the motion of the node and the law of change of inclination. 
He shewed that when Sun and Node were in conjunction, 
then for nearly a month the Moon moved in a plane very 
approximately, and that the inclination of the orbit then 
reached its maximum, namely, 5° 17' about; but as the Sun 
moved away from the Node the latter also began to move, 
attaining its greatest rate when the separation was a quadrant, 
and that at this instant the inclination was 5° very nearly. 
He also assigned the law for intermediate positions. The 
fact that there was no motion when the Sun was at the 
Node, that is, in the plane of the Moon's orbit, confirmed 
his theory that these inequalities were due to the Sun's action; 

When we spoke of Newton's results as fragmentary and 
incomplete, let it not be understood that he gave only very 
rude approximations to the truth. His results are far more 
accurate than those arrived at in elementary works of the 
present day upon the subject. 

After Newton, Clairaut (1713 — 1765) treated the Lunar 
Theory analytically. He readily found the Variation and many 
other inequalities, but met with a difficulty in determining the 
motion of the apse. At first he made its mean motion only 
about one-half of the observed value, and supposed that this 
indicated a failure of Newton's law of the inverse square of 
the distance; but soon he recognized an error, caused by 
omission of terms which he had imagined would not aflfect 
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the result. When these were included the calculated amount 
was nearly doubled. 

The first Tables of the Moon which were sufficiently accu- 
rate for use in determining longitudes at sea by observation 
of Lunar Distances were those of Mayer. They obtained a 
prize oflfered by our Board of Longitude, and were published 
in 1770 by Maskelyne, the Astronomer Royal. 

The first Theories which gave the Moon's place with an 
accuracy equal to that of observation were those of Damoiseau 
and Plana. The former was published in 1827, preceded in 
1824 by Tables; the latter was published in 1832. 

Hansen's Tables, which are those now used, were constructed 
from theory and were published in 1857 at the expense of the 
British Government.] 



LECTUEE II. 

ACCELERATIONS OF THE MOON RELATIVE TO THE 

EARTH. 

When three bodies move under their mutual attraction, 
their motions are unknown to us except in the cases when 
they are approximately elliptical ; but this restriction includes 
almost all the most important cases in the Solar System. 

If one body of the system is greatly predominant and if 
the lesser bodies are not close together, the centre of gravity 
of the greater body may be taken as a common focus around 
which the others move in approximate ellipses. Or again, if 
two bodies lie close together, their relative motion may be 
approximately the same as though they were isolated, although 
the system contains a third greatly predominant body; for 
their relative motion is aflfected by the difference of the attrac- 
tions of the central body upon them and not by the absolute 
value of those attractions. 

The Sun and Planets are an example of the first kind; 
the Earth, Moon and Sun of the second. The Earth and 
Moon describe orbits round the Sun which are approximately 
ellipses, and the Moon might be regarded as one of the. planets; 
but this point of view would not be a simple one ; the disturb- 
ing action of the Earth would be too great, though it is never 
so great as the direct attraction of the Sun, that is to say, 
never great enough to make the Moon's path convex to the 
Sun. The more convenient method is to refer the motion of 
the Moon to the Earth, and counting only the diflference of 
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the attractions of the Sun upon the Earth and upon the Moon, 
to find how this distorts the otherwise elliptical relative orbit. 
This is the method of the Lunar Theory. 

The position of the Sun must be referred to the same 
origin ; but since the Earth describes an ellipse about the Sun 
which is disturbed by the action of the Moon> if we choose as 
origin the Earth's centre, we must allow for the disturbance of 
the Sun's position by the Moon. This correction may be evaded 
by choosing as origin, not the Earth's centre, but the centre of 
gravity of the Earth and Moon, 
with respect to which the Sun 
describes a curve so closely el- 
liptical that no allowance is 
required. For, if fif, E^ M denote 
respectively the Sun, Earth, 
and Moon, and the centre 
of gravity of E and M, the 
accelerating forces of S are 

on E 8/8E^ in E8, 

on Jf 8/8MHnMS; 
and these imply accelerations of of amount 

E 8 

^^:jr^ ^^. parallel ^ MS; 

now the accelerations of S are 

EjSE^ in SE, 
M/SM^ in SM; 
hence the acceleration of relative to /S is 

8 + E + M E „ 1 + po 

S-\-E + M M 1, 1 , M-« 

-^tm- sm P*"*"^^ ^ ^^' 

S + E+M/„ OE „ GM\ . ^„ 
E + M r-SE*-^-8M») '"^ ^^' 

S + E + M („ SO , ^ S0\ . ^„ 
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Let EM=^r, 8G = r\ SOM^^to] then 

E M 
QM- r 0E=—— r 

Hence 



1 1 r jCi T 

]^, = /i [1 + -ETT-lr - 3 cos 6) 



E r 



S 



, / E rVf 3 . 15 , \ . "I 

1 1 r, M r „ 

77t^ =^ ~r^\ ^'- -ri i> —/ O COS 0) 

SE^ /^L E + Mr 

( M rVf 3 15 , \ "I 

and the accelerations of G are 
S + E + MV EM 7^^ ^ ^■ niur 

Now r// is approximately j^ ; neglecting the square of this 
quantity, we see that S moves about (r in a pure ellipse. 

Consider now the accelerations of the Moon relative to the 
Earth ; subtracting the accelerations of the Earth from those of 
the Moon, we find 

E^M ^/MO . E0\ . ,,^ 

let E'\-M=fi, S=^m ; then these become 
/A , mVf- , E-M r ^ , "I . 



inJf^, 



tVf^ ^E^Mrf 3^15 , V 1 

7a[3cos(. + ^^^^^-2 + yCOS^(.j + ......J 



parallel to OS, 
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[lect. 




In the accompanying spherical 
triangle, let be the centre of the 
sphere, 8M' the ecliptic, and M' 
the projection of M, 

Let 1/u be the projection of £ 
ME on the plane of the ecliptic ; 

the longitude of the Moon as seen from the Earth, 

(y the longitude of the Sun as seen from 0, 

s the tangent of the Moon's latitude MM\ 

Then 

r = (1 + sPf w\ COS G) = cos {0 - ff) (1 + s^)-^, 
and the accelerations of M relative to E are 



m 



+ — ^^ — - 1 + ET-r-Tf D- 3 cos (d - ^ + m Jf£^, 

1+52 r'ia L E + Mr^u \ ^^ 1 



/ii^'' 



,]. 



m 



'»« |_ ^ ^^ E + Mru 



parallel to OS. 



Call these quantities U and F respectively ; then if we 
resolve parallel to Jf' (?, perpendicular to M'O in the plane of 
the ecliptic, and perpendicular to the plane of the ecliptic, we 
have the following quantities which we call P, T, S; viz. : — 

P = U{1+ s')-^ - Fcos (0 - 0'), 
r=-Fsin(d-d')> 

and also 

Sf-P5= Fscos(^-^'). 

From these we find 

-, fiu' m' fl 3 ^.. ^, . E^M 1 



(^g-|«^)cos(d-eO + ^cos3(^-5') 



+ ... 
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/«w L 2 ^ ^ E+Mru 

|(|-|5^)siii(d-^) + ^sin3(d-^)| + ...], 

|(f-|«')cO8(d-^) + ^C083(d-^')} + -]- 

Hence with the time as independent variable we have the 
equations of motion 

d^ \dt) ~ ^' 
rdt\ dt) 

Or we may write these with as independent variable ; let 

SO that -jT = iiu\ 

dt 

Then H -j^ = — , 

ad li' 

whence ^*«'(d^ + «j + ^ W« 5^ = "^' 



agam, 



whence ^•.«.(g + «) +^f .^-P.-S; 
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or the equations of motion may be written 



H-«-l^..) = F-T^. 



TrdH_T 
d0 ~v?' 

^■»'(S+«)=-p'-«-^S- 

Our problem is to discuss these equations and to obtain 
from them expressions for the Moon's position at any time. 
The integration is best eflTected by observing what kinds of 
terms will disappear on substitution in the equations, and then 
assuming for the desired expressions for the coordinates a series 
of such terms multiplied by undetermined coeflScients. Our 
procedure will be to discuss one by one the irregularities which 
can be isolated from one another. This will permit a survey of 
the entire field without involving needless complexity; but if 
the Lunar Theory is to be accurate, the combinations of such 
terms with one another must also be included, and the number 
of terms employed and the labour of manipulating them becomes 
very great. 



LECTUEE III. 



THE sun's coordinates IN TERMS OF THE TIME. 



To obtain the Moon's coordinates in terms of the time from 
the equations found in Lecture II., we must substitute in the 
expressions for the forces the developments of the Sun's co- 
ordinates which we now proceed to give. 

Employing as coordinates /, ff, of the last lecture, we have 
seen that the Sun's motion may be regarded as purely elliptical, 
so that 

g' _ 1 + e' cos {& - istQ 

^ - tsr' = ri'« - isr' + 2e' sin (n'e - isr') + ^ e'« sin 2 (w'« - tsrO + . .. 

in which we have written for convenience vSt in place of v!t + e'. 
The quantities that enter the equations are 

/a' V cos g. .^ >wv 
\r/'sm ^ ^ 

Making the substitutions we find without difficulty 
(^Y= 1 + |e'' + 3e'cos(w'«- w' + |c'«cos 2(n'«-«r') + ... , 
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\r / Sin ^ ^ V 2 / sm ^ ■ 



+ 

,'\4 



\r / sm ^ ' ^ '^ sm ^ "^ 



+ 3«' T {(^ - n'O - (n'« - isr')} 
sm ^^ ^ ^ ^^ 

sin ^^ ^ ^ ^ 

8 sm i^ ^ ^ ^> 

+ y e'« ""^^ {(^ - n't) + 2 (n'^ - tsr')} 
8 sm i^ ^ ^ ^^ 



f ^' V ''!'^ 3 (d - ^) = (1 - 6e'0 """^^ 3 (d - n'O 
\r / sm ^ ^ ^ ^sm ^ ^ 

+ Sc' "^P^ {3 (6 - Ti'O - (^'* - ^0} 
sm ^ ^ ^ 

sm ^ ^ / V /J 

+ l|Ie'«T{3(d-n'0-2(V«-tsrO} 
8 sm V 

+ i 6'2 ^^^ {3 (^ - n't) + 2 (7i'^ - ^')\ 
8 sm ^ ^ ^ ' 

+ 

These quantities are to be substituted where they occur in 
the expressions for the forces found in Lecture II. 

Let us now make a few general remarks upon the result of 
the substitution. 
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It will be observed that the disturbing forces all involve the 
coefficient mV""'. It is very important to notice that the Sun's 
parallax is not required for the evaluation of this quantity. By 
Kepler's Third Law it is derivable from observations of the 
Sun's mean motion alone. Other terms however, namely those 
with the coefficient rn/a'Ht, involve the Sun's parallax directly ; 
and that constant may be obtained by comparing the observed 
with the theoretical values of the coefficients of those in- 
equalities, with an accuracy probably greater than that of any 
other method. 

The mean disturbing force is radial, and is equal to 

or the mean eflfect of the Sun's disturbance is to diminish the 
Moon's gravity towards the Earth ; and to diminish it more, the 
greater is the eccentricity of the Sun's orbit. Now e' has been 
diminishing for ages; hence the Moon's gravity towards the 
Earth has been increasing, and its average time for accom- 
plishing a revolution about the Earth has been diminishing. 

This is one cause of the Secular Acceleration of the Moon's 
mean motion which Halley derived from the records of ancient 
eclipses. 

It may also be noticed that the coefficient of the chief 

periodic part of the disturbing force, which involves 1 — „ ^'^ 
increases as e diminishes. 

Finally let it be observed that the term with argument 

2 (^ - n't) + 2 {n't - tsr'), 
which does not involve the Sun's Mean Longitude, is absent 

from the development of f— j . 2(0 — ff). 



LECTURE IV. 



THE VARIATION. 



The Variation is the first inequality we shall consider ; this 
is the inequality which is independent of eccentricities and 
mutual inclination in the orbits of the Sun and Moon. 

Let us first take the equations in the first form in which 
they are given in Lecture II., namely with t as independent 
variable : 

r dt V dt)' ' 

we omit the equation of motion in latitude, and in the ex- 
pressions for P, T we suppose « = 0; moreover it is possible 
and convenient to discuss separately the terms that involve 
the Sun's parallax; let these be omitted and we have 

rd^- U)+H=2/-. + 2/i«°«2<^-^>• 
and if c' = 0, r^ = a', m'/a'^^n'', 0' = n't + e, 

dt^ + rdt^ = -2n'8in2(^-n*-e); 
these are the equations to discuss. 
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Assume as a first approximation 

= nt-\'€ + b2ain[2(nt + €)-2 (n't + e)} 

= 71^ + e + &a sin 2-^, say ; 

1 1 

- = - [1 + Oa cos 2Vr], 

and we shall suppose Oa, b^ so small that in the first instance we 
may neglect their squares and products. 

Substitute in the equations ; then 
4 (n —ny ag cos 2-^ — {n^ + 4n (n — n') b^ cos 2^} 

+ -3(1 + SOaCOS 2i/r} = 1 n'^ + 1 n'»cos 2i|r 

3 

— 4 (n — n )2 &2 sin 2*^ + 4 (n — ri') r^Oj sin 2-^ = — s ^'^ sin 2-^. 



Hence, equating the coefficients of similar terms, we have 

'i = n' + ln", 



a« 



which gives the relation between n the Moon's mean motion, 
and -, the mean of the reciprocal of the distance; also 

r4(n-n')« + ^"|a,-4w(n-w')62= |ri'' (1), 

- 4 (n - n^b^ + 4n(n - w')a2 = - 1 w « (2). 

From (2) 4n (n — w') 60 — ^f^^H = « / . 

Add to (1), and substitute for /t/a'; 

3 „ 2w-n' 1 

02 = s w • 






A. L. 



*t=* 






3 ^ « ii-i' I 3 1 



«^:= 1'- 



Ouiznq — =: !», w^ laTii 

3 . 2-» 1 



Hi 

, 3 , 2-» 1 3 »* 



± ll-jnr" 11 . >il-llir^ 

3 — '5« — -^MT 



#ir, GkTlia^ , , or = i». , we !iii7e 



» — »" 1 — III 



% = 5"*l'- 



^--l - 1 

3_, a-hw,>'2-nii> . »_. 
&.=2"*^- i — '^s*'- 

3-2«,-hi|iH* 

TWeat aie eoiiTenient expresnoas^ ausd, as it happens^ Teiy 
ayproTTmaie, If we wish to derelope in ascendiiig powers of 
m or w^i, it ap|)€ais that the lacter deTelopcnent will be the 
DMte coQTeigent. 

We find hy obBeiTati*~>o — = "0748*), Tety nearly. 

Hence a, =: ^)0T17 J95, 

i^ = 01021^ = 210«"-4. 

Hence the ratio ^ the greatest and least distances will be 

1-00717^5 : 099282,05, 

and the greatest angnlar deviation from the mean IcMigitade 
win be 

3.5'e"-4. 

a "s^ dose approximation to the truth. 
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Also we have found 

= 71^ X 1-00280, 

which is the relation between the actual mean motion and the 
actual mean distance (or rather mean reciprocal distance) of the 
Moon. 

Without the Sun's disturbing action, the relation between 
the mean distance and the mean motion, or rather between the 
radius of the orbit supposed circular and the uniform rate of 
angular motion along it would be 

a* 

Hence in the actual orbit, the mean motion for a given 
mean distance is smaller than it would be without disturbance ; 

Or, for a given mean motion, the mean distance is smaller 
than it would be without disturbance. 

In fact, the relation between the mean distance and the 
mean motion is the same as it would be if the sum of the 
masses of the Earth and Moon were diminished in the ratio of 
1-00280 to 1. 



2—2 



LECTURE V. 

THE VARIATION {continued). 

I 

We will now proceed to substitute in the diflTerential 
equations the values of 1/r and which we have obtained, 
retaining terms of the order of the squares and products of 
Oa, &2 and m' or mj^ 

The values to be thus substituted are 

- == - (1 + Oa cos 2^), 
r a^ ^ 

^ = n^ + 6 + &a sin 2-^, 

where yjr^nt + e — (n't + e'), 

3 , 2 + m, 



Hence 
r 



rdt^ 



3 

62 = (1 + ^) <h + g Wi*. 



= a 1 — a2Cos2'^ + ^aa^(l + cos 4^) 
= 4a (n — ny [a^ cos 2^ -- 2a2^ cos 4-^], 

ri 3 

= 4 (n — ny \a^2^ + (h cos 2^ — ^ ag^ cos 4-^ 
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again, 



-IT = w + 2 (n - w') 6a cos 2i/r, 



dt 
[-it] = n^ + 4w (n - n') 62 cos 2i^ + 2 (n - nyh^ [1 + cos 4^], 

^= Jl + So^cos 2f + |aaHl + cos4np^)J . 
Also, 
- ^ = 2 (n - w') Oa sin 2^ - 2 Oa^sin*^ , 

- -TT -^ = 2(n— n') noa sin 2-^ + ](n — n') Oaia — s ^^' • sii^ ^^^ > 

-^ = — 4 (n — w')' 6a sin 2-1^. 

And cos 2 (^ — n'^ — e') = cos 2^ — 6a (1 — cos 4-^), 
sin 2 (^ — 71'^ — eO = sin 2-1^ + 6a sin 4-^. 

Substitute these in the difiFerential equations, and we get, 
on transposing all the terms to the left-hand sides from the first 
equation 

^{n — ny ^02^ + OaCos 2'i^-^a2^cos4'^ 

-[nH2(n-n)"6aH4n(7i-n')62Cos2'f+2(n-n76a'cos4Vr] 

+ ^3 1 + 2 02'+ SOa cos 2^ + 2 Oa^ COS 4^ 

1 3 

- 2 n'^ - 2 n'2 [- 63 + COS 2'^ + 6a cos 4i^], 

and from the second equation 

- 4 (n - ny 62 sin 2'^ 

4 (n — n') wog sin 2-^ + \{^ — '^0 ^2 — 5 ^cta'f sii^ ^^^ 



+ 



3 
+ 5 w'" [sin 2-1^ + 62 sin 4i/r]. 
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The coefiScient of cos 2-^ in the first of these expressions, 
and that of sin 2i^ in the second, are respectively 

4(n - nja^- in (n -ti') &2 + 3 ^ a, - ? n\ 

3 

and - 4 (n - rif h^ + 4n (n - v!) a, + ^ n'«, 

and these are evidently reduced to zero by giving Oa, \ the 

values previously found, if we substitute for /^t/a' the approxi- 

1 
mate value w* + « ^^'"* To find the more correct value of /tt/a', 

equate to zero the constant term in the first expression ; 

2(n-n')»a.»-n«-2(w-n W+g (l + |a2^) -^ n'»+ |n'>6, = 0, 
that is 

1 
Hence we see that /tt/a' differs from w^ + ^ v!^ only in terms 

of the fourth order, if we consider rn^ a quantity of the first 
order and consequently Og, 62 quantities of the second order. 
Hence also by taking 

5 = n' + i«^ 

in the multiplier of a^, when we equate to zero the coefiScient 
of cos 2^, we only negle(fted a quantity of the sixth order in mi, 
and the error in the resulting values of Oj, 63 is of that order. 

We see that the substitution just made in our equations 
leaves outstanding terms of the fourth order in cos 4i^ and 
sin 4^. In order to get rid of these we must add terms of this 
form to the assumed values of 1/r and B, respectively. Suppose 
that 

- = - [1 4- Og cos 2^ + a4 cos 4-^], 
^ as n^ + € + 6a sin 2^ + 64 sin 4i^, 
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where^ as we shall find, a4 and 64 are small quantities of the 
fourth order. 

It may be readily seen that the additional terms introduced 
are the following : — 

1 d?r 
in -^ 16(w-n7a4Cos4^, 

— I -J- j — 8n (n — n') 64 cos 4-^, 
^ ^ 3a4 cos 4-1^, 

-^ - 16 (n - nf 64 sin 4^, 

rdtTt 8«(n-n)a«sin4t. 

and also that we may neglect the terms added in the ex- 
pressions for 

|n'«cos2(^-n'e-€'), | n'« sin 2 (^ - n'^ - e'). 

If we write the terms thus produced along with the several 
terms lefb outstanding in our equations and then equate the 
whole to zero, we have 

16(n-n')'a4-8n(M-n')&4 + ^3a4-6(n-nO»a,^-2(n-n')»6a« 

-16(n-ri764+8w(n-n')a4+4(ri-n7aA-2n(w-7i')aa" 

bom which we must determine a^ and 64. 

Put -8 = ^' + ^^^ 

a' 2 

g 
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and divide both equations by (n — n')'; we get 

3 / 3 \n 27 

-2(^l+2mi+2^2jla,»-3(H-ma)mi»aa-32^i*=0, 

8(H-mi)a4-1664 + [4(l+mi)-2(l+mi)]a22 

(1 \ 9 

Simplify and multiply the last equation by ^ (1 + mi), 

(9 \ /13 1 \ 

19 + 6mi + 2 7Wi«ja4-8(l + mi)64-(-2+mi--|mi«ja22 

27 
- 3 ( 1 + mi) mi« Oa - g2 W = 0, 

(4 + Smi + 4mi2) a* - 8 (1 + mi) 64 + (1 + 2mi + mi^) a^ 

3 / 1 \ 9 

+ 2(1+^1) (^1 +2mijmi2aa + g^(l + mi)wii*=0. 

Subtract the latter from the former and 64 will be elimi- 
nated ; we get 

^15 - 2mi-|-2 mi3Ja4 " (y "*" ^^ "*" 4 ^'') ^' 

3 / 1 \ 9 

-g(l + mi)f3 +2mijmi3a2-g2(4 + ^)^i* = 0, 

which gives a4 ; and this being known 64 is found from 

64 = 2 (1 + Wi)a4+ 3 (1 + rn^)(^% "^ leV^ "*" 2 ^) ^'"^ "^ 256 ^** 

Taking m = '07480 as in Lecture IV, we find 

a^ = -00004,580, 
64 = 00004,237 = 8"-740. 



LECTURE VI. 

THE TABIATION {continued). 

Let us consider the problem of the Variation over again, 
taking now as independent variable. 

The equations of motion are given in Lecture II : — 

0^ _^P T_^ 

dH T 



H 



dd u* 



where _ =^__ _ -_ _ cos 2(^-^), 

^= -|^8ia2(^-n 

u^ 2 It* ' 

80 that the second equation may be written 

Our aim is to express t and u in terms of and constant 
quantities. Now since the orbit of the Moon does not differ 
widely from a circle we may write the difference of nt + e from 
6y and the difference of au from unity as series of small periodic 
terms depending upon 0. Inspecting the form of the equations, 
it is evident that these periodic terms are of argument 2 (^ — ^) 
and its multiples ; that is 

7i^ + e = ^ -f periodic terms of argument 2 (^ — ff\ &c.; 

but n't + e' = ^' ; 
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therefore 

^-^ = (l-m)^-)8 

+ periodic terms of argument 2 (1 — m) ^ — 2/8, &c., 

where we have written 

^ = e' — wi€ ; 

this constant /8 is associated with (l — m)0 wherever the latter 
occurs ; for brevity in writing, we shall omit it. 

We may then assume as a first approximation 

aw = 1 + Oj cos (2 — 2m) 0, 

Tie + € = ^ + 63 sin (2 - 2m) 0] 
whence 

2(0-0')^(2-2m)0- 2mb^ sin (2 - 2m) 0, 
cos 2 (^ — ^) = mJa + cos (2 — 2m) - mh^ cos (4 — 4m) 0, 
sin 2 (^ - ^) = sin (2 - 2m) ^ - m\ sin (4 - 4m) 0, 

II ^ = 1 + (2 - 2m) 62 cos (2 - 2m) 0. 

Substitute in the right-hand member of the second equa- 
tion : — 

^^ ^?^|^ = - 3m« [sin (2 - 2m) ^ + (2 - 3m) 62 sin (4 - 4m) e]. 
Therefore 

which we may write 

log^ (C') = 2A2 cos (2 -2m) 6 + 2K cos (4 - 4m) 0, 

where A is an arbitrary constant of integration, A2 is a known 
quantity, and h^ involves 6a. If we take as a second approxi- 
mation 

aw = 1 -h Og cos (2 — 2m) + a4, cos (4 — 4m) 0, 

71^ + e = ^ + 62 sin (2 - 2m) ^ + 64 sin (4 - 4m) 0, 

the above value of logeiH^fh^) will not require modification 
and will supply equations of condition for determining the 
coefficients Oa, 62, a^, 64. 
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rr^^ dt U Utt^ k 1 

Thus yi -77i = 



d0 Hv? h H(auy* 
so that 

but 
loge (n ^) = -(1 -m)»6,» + (2- 2m)6,co8(2 - 2m) 

+ [(4 - 4m) 64 - (1 - my 6j»] cos (4 - 4m) 0, 

loge at* = — -J- -f Oa COS (2 — 2m) ^ + f a4— -|^ j cos (4 — 4wi) ft 
Hence we find 

(2-2m)62 = -A2-2a2, 

1 
(4 - 4m) 64 - ( 1 - m)» h^^ = - A4 - 2a4 + ^ a,'. 

The remaining equations of condition that we require are 
obtained from the first equation of motion ; this may be 
written 



d^ (au) 



hi("'S'-i(»'S)"-^<^-H 



d(ar.)3^^,|J^^2(^-ft) = g. 



dd 2 
Now 

aw = 1 + Oa cos (2 — 2m) 0-\' a^ cos (4 — 4m) 0, 

whence 

-^^ = - (2 - 2m) Oa sin (2 - 2m) ^ - (4 - 4m) a4 sin (4 - 4m) 0, 

^^ = - (2 - 2m)X cos (2 - 2m) ^ - (4 - ^ifa, cos (4 - 4m) 0, 

and 

(n'S =m''[H-(4-4m)6,C08(2-2m)^], 

n' ^Ysin 2{0-e')= m" [ain (2-2m)^+(2 - 3m)6j 8in(4-4m)^j, 
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= m« [(2 - to) &, + cos (2 - 2m) ^ + (2 - 3m) b^ cos (4 - 4im)0], 
^, = ^[l+V-2AjCos(2-2m)^ + (V-2A,)co8(4-4m)^. 

Substitute these in the equation above, and equate the 
coefficients of corresponding terms, 

l + |m» + |m''(2-m)6,+ |m^ + |TO'(2-2TO)a, = ^(l+V) 

- (2 - 2m)«a, + (l + 1 TO«) a, + 1 m» + (2 - 2m) m«6, = ^ (- 2A.) 

(1 \ 3 3 

1 4-2mMa4 + 2w'(2-3m)6a+ jm^Oa 

If we neglect at first terms of the fourth order, we find fi*om 
the first of these equations 

From the earlier set of equations we have 

(2-2m)fe2 = -Aa-2a3; 
substitute this in the second equation above. We get 

-(2- 2m)«+ l+|m2- 27/42^ - w»Aa + |m2 = (l + 2m«)(--.2A,), 

or 

f3-8m + M^')a. = |m' + 2A.. = |m' + ?-^-- = ?m«^::^. 
\ 2/2 2 21— m 2 1 — m 

so that 

_3 , 2-m 1 

'^""2'^ r^^ ^ . 11 ' 

3-8m + ^m2 



and 



^- 1 « 1 t 

'^'"~l-m*^"2^^2^'^ 
1 Z m^ 



1-m^ 8(l-m)''" 
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These are numerically equal to the quantities denoted by 
the same symbols in Lecture IV, but 6a bears the contrary sign. 
We further find 



(l5-32m + YmAa4 



= 7 r^ {(1 + 3m - 2m») a, + (8 ^ 15m + 6m«) b,l 
4 1 — m ^ ' « \ /J 

or reduced to numbers 

a4 = - 00002,210, 

64= -00005,414 = in7. 

Finally let us exhibit the relation between the constants 
employed in this investigation and those of Lectures IV, V ; to 
distinguish them, attach accents to the latter, so that 

^ = 71^ + 6 + 6/ sin 2-^ + 6/ sin 4^, 

a' 

— = 1 + da' cos 2i^ + a/ cos 4^, 

and, omitting the constant /8 as before, 

(1 - m) ^ = 1^ + (1 - m) bz sin 2i^ + (1 - m) b/ sin ^yfr. 
Then 

2^ = (2- 2m)^- (2 - 2m)6;sin(2 - 2m) 0, 
sin 2^ = sin(2 - 2m) ^ -(1 - m)^/ sin (4 - 4m) 0, 
cos2^ = (l-wi)6; + cos(2-2m)^-(l-m)6/cos(4-4m)^. 
Substitute in the equation for ; we find 

7i^ + € = ^-62'sin(2-2m)^-[64'-(l-m)6a'2]8in(4-4m)^, 
and similarly 

a'l^ = 1 + (1 - m) Oa'Ja' + aV cos (2 - 2m) 

+ [a4 — (1 — m) Oa'fta'] cos (4 — 4m) 0. 

We observe that a' diflers from a by quantities of the fourth 
order. 



LECTURE VII. 

CORRECTION OF APPROXIMATE SOLUTIONS. 

We may simplify the equations we have been dealing with, 
by a proper choice of units. Let the unit of distance be the 
radius of the circular orbit which the Moon, if undisturbed, 
would describe about the Earth in its actual periodic time; 
then 

Also choose the unit of time so that 

n — n' = 1, 

so that, if we take as the result of observation of the mean 
motions of the Sun and Moon, 

n '.n^ -07480,13, 
we get n' = -08084,9 = mi, 

where rrii is the quantity so called in Lecture IV; and 

/^ = 1-16823,4. 
We shall frequently adopt these simplifications in what follows. 

Now let I = loge {rla)y 

so that 

and the equations discussed in Lecture IV become 

dH (dly (d0\^ fi ^ ,. fl 3 o "1 ^ 



^ dld0 

dt^ "^ dt dt 



n'^\ |sin2o>1=0. 



where a} = — 0\ 
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Now these equations are defective, for they have been 
formed by omitting certain terms from the complete equations 
as given in Lecture II. Hence, calling Iq, 0q the values of I, 0, 
which we have proved in Lecture IV to be solutions of the 
above equations, if we substitute l^, 0o in the complete equa- 
tions of Lecture II, residuals are left, say X and Y respectively. 
And if I, be solutions of the complete equations, and if we 
write 

= 00 + 80, 

where 81, B0 are small quantities whose squares and pi*oducts 
may be neglected in the first instance, we obtain the following 
equations for determining 81, 80, the corrections to approximate 
solutions lo, 00 already found : — 

Now let us write 

where c = (1 + nj + ^ n\ 

The quantity v consists wholly of periodic terms of the form 
cos 2iyfr multiplied by small coefficients ; w contains, besides 
periodic terms, a small constant term, which however might be 
removed if we were to choose c as the constant part of fi/vo^ in 
place of according to the definition above. 

Let 8^1, 8'0 be quantities defined by the equations 

then 8% 8'0 are approximations to the complete corrections 81, 
80, which if substituted in the equations that give those cor- 
rections will leave residuals, say X' and Y\ where 
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^' = 2 §« ^ - 2i; ^ - SwS'l + 3n'» sin 2o,S'^, 
at at at 

T' = 2^^+ 2v ^ + 3n'» cos 2a,SU 

We see that their value is known when S7, S'0 are determined. 

Now S7, B'0 may be determined as follows. 

Let X =Pq + 2pi cos i^, Y = Sg'i sin i-^, 

where i takes all positive integral values ; and assume 

S7 = ao + Sttf cos i^, S'^ = 26»- sin i-^. 

Then substituting and equating coefficients, the constant 
term gives 

Pq - Scoo = 0, 

and the terms in i*^ give 

Pi - tX- - 2 (1 + nO i6i - Scdi = 0, 
qi-i%-2(l+nf)iai =0; 

the second of these may be written 

2(l + r?')*-4(l+w7ai-^2(H-n')i6i = 0; 
subtract from the first and we have 

V 

^,-2(l+«')f 

or Oi — ———^^— 



and 6i = -2(l + n')"* + |. 

We see that «»♦, 6< will be of the same order of small quantities 
as jpi, g'i, in general. And therefore the coefficients of the terms 
of X\ T* will be of order higher than those of X, F. Proceed 
then to determine further corrections S"Z, i"0 satisfying the 
equations 
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then if S'Z + S'7, S'^ + O are substituted in the complete 
equations for Bl, SO the residuals become 

X" = 2 ^ ^^ - 2t> -^'^ - 3wS"l + 3n'» sin 2a,0, 

Y" = 2§^ + 2v^^'/ +Sn'*coB2<oB"0. 

at at at 

expressions which, if developed in series of cosines and sines of 
multiples of ^, will have coefficients of higher order than the 
corresponding coefficients in X\ Y\ The like process may be 
repeated until the residuals become insensible; we then have 
sensibly connect values of SI, S0, giving 

i = Z,+gi = Z, + S7 + S'7 + , 

^ = uq + ou = ^o4'5ff4"S -{- • 

We may now take into account squares and products of the 
small quantities SI, SO by treating lo'\-Sl, 0q'\-S0 as given 
approximate solutions just as we have here treated Zo> ^o; 
substitute them in the complete equations of motion, and 
determine the residuals X, Y which they leave. These resi- 
duals will form the basis of a second approximation, and the 
operation may be repeated until no further correction is 
necessary. It is to be observed that if SI, SO depend upon 
some such constant as the eccentricity of the Earth's orbit 
around the Sun, or the parallax of the Sun, then successive 
approximations yield correctly and separately the terms which 
depend upon the first, second, ...... powers of that constant. 



A. L. 



LECTURE VIIL 

THE PARALLACTIC INEQUALITY. 

We shall now apply the method of the last lecture to find 
the terms in the Moon's coordinates which depend upon the 
parallax of the Sun. 

The values of I, found in Lecture IV are 

Zo = \oge (r/a) = - Oa cos 2^, 

0o = nt -\- e+b^ sin 2^, 

and these satisfy the equations of motion in which the terms 
involving the Sun's parallax are omitted. Hence the residuals 
they leave from the complete equations are 

Z = -X7i'«^|^cos(^o-^) + ^cos3(^o-^')j, 



F= Xn'«-|i sin 
a 8 



(^o-^') + y8in3(^o-^)}, 



, ^ E-Ma 

where X= ^i — i>-7« 

Now from above 

^o-^' = '^ + 63 8in2^; 
hence we have 

sin (^0 - ^) = sin -i/r + - 62 (sin -i/r + sin 3^), 

COS(^o-^')= cos 1^ - - 62 (cos -l/r - cos 3'^), 

3 

sin 3 (^0 - ^) = sin 3'^ + ^ 62 (sin '^ + sin S-i/r), 

cos 3 (^0 - ^) = cos 3^ - ^ 6a (cos '^ - COS 5'^); 



+ 
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and 

/15 9 , 9 \ ^ , /45 , 15 \ 
l"8'"*'i6 ' " Te'^V^^ ''"*■ V16 ^«- ie'^j^^^^'''' 

v¥ 16 '^ " l6 V ®^ '^ V16 ' "^ 16 7®^^ ''''^* 
Assume 

— S\ = \ai cos ^ + Xoj cos 3i/r, 

S^ =: X6i sin -^ + \bi sin 3'i/r, 

neglecting for the present the tenns in 5-^. In the present 
case it happens that it is more advantageous to substitute 
these expressions directly in the complete equations for Bl, S0 
given in the last lecture than to follow exactly the process for 
finding them by successive approximation. Omitting the 
factor X, we get 

Oi cos yfr + 9aj cos 3-^ + 4a2 sin 2-^ [aj sin yfr + Sa^ sin 3-^] 

+3— 3[aiCos'^+a8Cos3i/r]+3^[3aaCos2'^][aiCosi^+a8Cos3'^] 

-[2(1+ n') + 462 cos 2i^] [61 cos yfr + 36, cos 3-^ ] 
-f 3w'* [sin 2-^ + 6a sin 4-^] [fh sin -^ + 63 sin 3-^] 

„ /9 27 , 3 \ , , /15 9 , 9 \ 

(45 15 \ 

jg 62 - jg ^2 1 cos 5'^ = 0, 

— 6i sin -^ — 963 sin 3i/r + 40, sin 2yjr [bi cos -^ + 36, cos 3^] 

+ 3w'2 [- 6a + cos 2-^ + 63 cos 4^] [6i sin ^ + 6, sin 3'^] 

+[2 (1 + n) + 46a cos 2'^] [oi sin '^ + 3a, sin 3'^] 

,,/3 21 , 3 \ . , ,,/15 3 , 3 \ . ^, 
V8 "" 8" ^ "" 4^/ ®^"^ "^ ^ l"8" ■*■ 16 ' " 16 7 ®^^ ^ 

'45 . 15 



■*" ^'' (iB *' " IB ^) ®^" ^^ "^ ^* 



3—2 



+ 
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If we equate to zero the coefficients of cos -^ and cos Syfr in 
the first, and those of sin yjt and sin 3^^ in the second, we obtain 
the following equations for Oi, 61, a,, 63 ; the terms in 5-^ remain 
outstanding, and the effect of 05, 65 in modifying the other co- 
efficient is neglected. 

a, [^1 + 2a, + 3 ^, (1 +|a,)] -6, |^2 (1 + n') + 26, - 1 n'^1 

a,[6a, + |^, J-6, r66,-|n"(l+6.)] = |n"(3-96,-4a,), 
a,[2(l + n')-26,]-6,ri + 1 «'» - 2o, + 3n''6,l 

+0, [66,] _ J, r6a,-|«'>+|n'>6,l = -|»'>(1 -76,-2a,), 

«! [- 2a, + 1 ^, a,] + 6. [- 26, - 1 »'» + 1 n'%,1 

+a,[9 + 3^,] _6.[6(l + »')] = |„'.(5+|6,-|a,), 

Oi[26,] + 6,r2a, + |n"-|n'«6,1 

+ a,[6(l+»')] -6,[9 + 3«"6J =-|'^''(5+|^>-|«»). 

If we require the formal values of Ui, 6i, a,, 6,, we must 
substitute for a,, 6„ /i/a' the expressions we have found for 
them, and it will then be best to develope the coefficients in 
ascending powers of n'. £ut it is difficult to obtain by this 
process such good numerical results as we can get by sub- 
stituting the numerical values of a,, 6„ fi/a' immediately in 
the equations above. If we do so we get the equations 

4o6672ai-2-172326,+ 080930,- 051376,= ?n'»x 2-87937, 

O 

214128ai-0-995646,+ 061270,- •033386,=-|»'»x 0-91416, 
-02349a,- •030126,+12-51451a,- 6*485086,= |n'»x500455, 
•020420,+ -024066,+ 6-48508a,-9000206,=-|n'»x 500152. 
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We notice that the first equation is not very different from 
the second doubled; it is this fact that makes successive 
approximation a disadvantageous method and renders it ad- 
visable to include small quantities from the beginning. 

Eliminate Os, 63 in succession from the third and fourth 
equations, thus : — 

Multiply the third equation by 

9-000204.[12-51451x 900020 6-48508 x 648508]= 0127523, 

and the fourth by 

-6-48508-T-[12ol451x900020-6-48508x6 48508]=-0091887, 

and add ; 63 will be eliminated. 

Again multiply the third equation by •091887 and the 
fourth by —•177317 and add; a^ will be eliminated. Hence 
we find 

•001 1 1 9ai - -00605261 + a, = | n'^ x 1-09776, 
- •001463ai - •0070346i + 6, = | n'» x 1-34669. 

o 

Multiply these by —•08093 and ^05137 respectively and add 
to the first equation : — 

4-56672ai-2-l72326i+08093a3-0513763=|n'»x 287937, 

o 

-000009ai+0-000496,--08093a, =|n'»x-0'08884, 

-0 00008a, -0-000366i +05137ft,=|n''x 0-06918; 

O 

hence 

4-56665ai - 2-l72196i = ^ n'^ x 285971. 

o 

Eliminate c/3, 63 in a similar manner frx)m the second 
equation ; 

2-14128ai-0-995466i+-06127a3--0333863=-|7i'«x 091416, 

o 

-000007a,+0000376,-06127a, =-tn''x 006726, 

O 

-000005ai-0000236i +-0333863=-|n'2 x -004495 ; 

o 
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hence 

214116ai-0-995506i=-|n''» x 093647. 
From these equations we find 

a, = - 1 n'2 X 46-4814 = - '11392,8, 

o 

6i = -?w''x 990336 =--24273,4, 
a,= |n"x -56042= -00134,9, 
6,= |m'»x -58209= -00142,7. 

O 



LECTURE IX. 

THE PARALLACTIC INEQUALITY (continued). 

Let us now consider the terms in dyfr which have been left 
outstanding. 

Include additional terms XOj cos 5yfr, X65 sin 5-^ in — Si, Stf, 
and equate to zero the coefficients of cos 5^, sin 5^ in the 
differential equations that give 81, 80. 

We have 
25a5 + -^ as - 10 (1 -h n')b6 - 60^^ - 6bJ><^ + ^ ^a^ 

- 2565 + 10 (1 + %') a, + 6a^, + 6a^, 

In these equations substitute 

o* 2 

Then 

(28 + 6»' + I «'») c, - 10 (1 + n') 6. = (I - 9n' - ^ n'») a^, 
-2ofc.flO(l+ii')a,=-6a^,-66^-|»'»6,-»'»(^6,-^«,). 
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Eliminate 65 : 
(24-2n' + ^n'«)a. = [^(|-9»'-^«'«)a, + ^(l + n')das 

+ I 66, + 2 n'> + ^1 (1 + n') a, + 1 (1 + n') n'A 6, 

and then 65 is given by 

From these we find 

a, = I n'» X -00595,3 = 'OOOOl 4591, 

6. = |n'» X 00710,3 = 00001,7410. 

O 

These numbers being so small, we see that we may safely 
ignore, as we have done, their efFect in modifying the earlier 
coefficients. 

To find the effect of these coefficients upon the Moon's 
coordinates we must multiply by the factor X = „ — =^ . — , 



We shall take in accordance with the results given in 
Monthly Notices, Vol. 13, p. 177, and Appendix to the Nautical 
Almanac, 1856, 

^ = 81-5. 
M 

Constant of Moon's Parallax = 3422"-325. 

Also we shall take in the first place, the Sun's Mean 
Parallax to be 8"-8, and in the next place 8"*9, and we will 
find the corresponding values of the coefficients of the Paral- 
lactic Inequalities. 



IX.] 
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We find 

8"-8 

X = -00250,9 
Xoj = - 00028,585 
X6, = - -00060,903 

Xo, = -00000,3385 
XAs = 00000,3580 

Xa, = 00000,00366 
XA, = 00000,00437 



\ = 



Xo, 
125"-62 X6i 

XOs 

0"-7384 X6, 

Xa, 
0" 00901 X6, 



8"-9 
•00253,76 

•00028,910 

•00061,596 
= -127"-05 

■00000,3423 

•00000,3620 

0"-7468 

•00000,00370 

•00000,00442 

0"00911. 



These results are very &irly accurate ; but in order to get 
good values for Oj, 6i, we were obliged to discuss o,, 6,, a,, 6, 
simultaneously. Let us consider the peculiarity of the equa- 
tions from which this difficulty arose. 

Following the method of approximation of Lecture VII, if 
we neglect at first the products of SI, S0, dBl/dt, dSd/dt with the 
small quantities Og, bj, n\ the equations become 



d^Bl _ ,^ dSe 
" dt 



dt* 
d^Be 



-2n"^.--3^8i + Z = 0, 
dSl 



a' 



dt^ ^^""dt 



+ F = 0. 



Now suppose the following is a set of terms that appear 
in X /)» cos (i^ + 7), in Y gisin(t^ + 7), 

SI ai cos (it + 7), B0 bi sin (it + 7) ; 

then as in Lecture VII, we find 



ai = 






6t = — 2-ai + -_g. 
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Therefore if i differs little from w, the divisor in Oi will be 
small, and a small error or omision in the numerator of di will 
appear magnified in the values of both Oi and 6{. In the case 
of the first term of the Parallactic Inequality, 

i = w — n\ 



and if we take 



i» - n» + o w '« = - 2nn' + ^ n'^ ; 



Pi = -gW'', ?t = gw'>, 



which differ from the correct values by quantities of the fourth 
order, then 

and the formulae give 

_3 ji'(5/i-3n') 
^* " 4 (?i -V) (4n - 5n') ' 

, _ 2n 3 n'« 

Now if we develope these expressions in ascending powers of m, 
i.e. n'/n, the first terms are 

15 , 15 

and these are the only terms which the formulae derived from 
our method of approximation will give correctly. 



LECTUKE X. 



THE ANNUAL EQUATION. 

Let us next take into account the effect of the first power 
of the eccentricity of the Earth's orbit. We shall find that it 
produces an inequality in the Moon's coordinates, the chief part 
of which has a period of one year, and is therefore called the 
Annual Equation. 

In the formulae of Lecture VII, let the known approximate 
solutions 2o, ^o» include the Variation only ; then the equations 
for the corrections 81, 80 are 



Z + ^J + 4a« sin 2Vr^- 3^,(1 + 3a-3 cos 2Vr)SZ 



- 2 [(1 + n') + 262 cos 2i/r] J^ + 3w'Hsin 2y^+ b^ sin 4^)8^ = 0, 
F-l-^ + 4a,sin2i^^ + 3ii'*(-62 + co8 2>;r + 63cos4^)Sd 
+ 2[(l+n') + 26,cos2Vr]^^ =0, 

where a,, 6a> /*/<** are known quantities whose values are given 
in Lectures IV, V. 

Refer now to Lecture III, and we find that the terms that 
are left outstanding when the terms of the Variation are sub- 
stituted, and the parallactic terms omitted are the following : — 
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+ I n V cos {2 (^ - n't) + {n't - w')} 

91 

F= + ^ nV sin {2 (tf - n'e) - {n't - tsr')} 

- I nV sin {2 (tf - n't) + (n'« - w')}- 

Write a for n'^ — tsr' ; then 
cos {2 (tf - n't) - a} = cos (2>;r - a) - (262 sin 2^) sin (2^^ - a) 

= - 62 cos a + cos (2-^ — a) + 6a cos (4-^ — a), 
sin j2 (5 — n't) — «} = + 62 sin a + sin (2-^ - a) + 62 sin (4^ — a), 
cos {2 {0 — n't) + a} = — 62 cos a + cos (2-^ + a) + 62 cos (4-^ + a), 
sin {2 {0 - n't) + a} = - 62 sin a + sin (2i^ + a) + 62 sin (4-^ + a). 

Hence 

X = -|7i'«(l-362)e'co8a-^nVcos(2i^-a) 

3 21 3 

+ -7 n'^e'co8 (2^ + «)- x w'^^^a^'cos (4>;r-a)+ - n'^/cos (4i^+a), 

F = 6n'*2c'sina + ^wVsin(2i^-a)-|nVsin(2>;r + a) 

21 3 

+ -J- n'^ftae' sin (4-^ " *) " 7 ^'^62^' sin (4^ + a). 

For our present purpose we shall ignore the small terms in 
4i^ — a and 4*^ + a which are of the sixth order. 

Assume 

— 81 = 056' cos a + a^' cos (2i/r — a) + Otc' cos (2'<|r + a), 

80 = 656' sin a + 6ge' sin (2-^ — a) + bje sin (2^ + a). 

Now the terms which arise in the left-hand members of the 
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equations owing to terms — Up cos pt in Bl, and bp sin pt in 80, 
will be 

p^Op cos pt + Za^pop [cos ( p^ — 2-^) — cos (pt + 2-^)] 

+ -^ ap\ cos pt + '^(h cos ( p^ — 2i^) + ^ Oa cos ( p^ + 2-^) 

- 2 (1 + n')pbp008pt - 2622)6p [cos (pt - 2>;r)+cos(|)* +2^^)] 
3 

+ 2 W''6p [cos {pt - 2'<|r) - cos (pt + 2>;r)], 

and 

— p^bp sin 2)^ + ^a^pbp [— sin (p^ - 2'<|r) + sin (pt + 2'<|r)] 

3 

+ o w'^ftp [sin (pt- 2^^) - 263 sin pt + sin (p^ + 2i^)] 

+ 2 (1 + »')pay sinp^ + 263^0^ [sin(2)^ -2'^)+sin(|)^+2^)], 
respectively, neglecting the very small quantities in 4-^. 

Hence we get the equations following : — 
Equate to zero the coefficients of cos a, sin a : 

(»'» + ^)a.-2/i'(l+n')6. + 2(2-»')a^ + ^.|a,a, 

-2(2-«')6A + |«'*ft. + 2(2 + n')o/i, + ^.|o^ 

- 2 (2 + n') bA + 1 n'% = | n'' (1 - 36,) 
- n'»6, - 3n'»6A + 2(1+ »') »'a, + 2 (2 - n') o^, 

- 2 (2 - n'),b^, - 1 «'«6, - 2 (2 + n') oA 

+ 2 (2 + n') ^7 + 1 n'% = - 6n'%' 
Equate the coefficients of cos (2^ — a), sin (2'^ — o) : — 

(2 - n')»a, + 3 ^, a, - 2 (1 + »') (2 - n') 6. + 2n'o8a. + 3 4 . 1 o«a. 

- 2n'hA + % n'\ = ^ »'». 



- (2 - n'f 6, - 3«"%A + 2 (1 + n') (2 - n) a, + 2n'aA - 2»'6^ 

3 /* 21 



- ? n'«6. = - ~ n\ 
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Equate the coefficients of cos (2^ + a), sin (2-^ + a) : — 
(2 + n)2a7 + 3 ^307-2(1 +n')(2 + n')67-2n'aaa5+ 3-3. la^as 

- 2n'bA - I n'% = - I n\ 

- (2 f n'y b, - 3»'»6A + 2 (1 + »') (2 +.n')a, + 2n'a,6, + 2n'5^ 



!«'*.= In' 



In equations of this class, as a general rule we would 
determine a,, 65 approximately from the first pair, substitute 
them in the second pair and determine a^, &« approximately, 
and similarly 07, b^ from the third pair, and repeat this approxi- 
mation as often as might be necessary. But if we refer to the 
second equation, we see that 65 must be determined by means 
of a small divisor, and this puts any method of approximation 
at a disadvantage! In order to obtain readily satisfactory values 
for the new coefficients, we shall treat the six equations simul- 
taneously, substituting first the numerical values of the known 
quantities. 

We have found 

as =00717,95, 62=01021,20, /i/a» = 117150,3. 
Hence 
3-52105a5 - 017476365 + •065405a« - -02939366 + 06772707 

- 03269567 = I n'» X 0969364, 

017476305 - OOO673665 - 03919706 + -01775366 + 04249907 

- 02007567 = I n'« X - 0040848, 

'039009O5 + 00815365 + 71976606 - 41485866 

= I n'« X 3-50, 

- 00165105- 00864365 + 41485806- 3-6833566 

= |n'*x-3-50, 
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•036687a,- •011455t, +7-8444.3ar- 4-498116, 

= |«'»x-0-50, 

•0016510,+ -0109656, +4-49811a,- 4-330126, 

= |n'»+ 0-50. 

From the second and third pairs we find 

•016186a, + -0070846, + a, = | «'» x 294739, 

•018678O, + -0103266, + 6, = | »'» x 4 26994, 
•011026a, - -0072036, + a, = | n'' x - 0-321398, 

•011072a, - -0100156, + 6, = | n" x - 0-449338. 

Eliminate a„ 6., a,, b, from the first pair. 
Hence 

3-52015a, - -1747626. = | w'» x 0-909170, 

0-1 748190, - -0065366, = I «'» x 0003516. 
Hence 

a, = I n'« X - 0-70619,8 = - -00692,37, 

6, = ? n'« X - 19-4268 = - 19046,3. 

Now e' is a constant found by observation; taking 
c' = 3459"-28, its value in 1850, we get 

o/ = - 00011,61, 

6/ = --00319,4 =-658"-9, 
and further 

a,= -03036,8, o^= -00050,9, 

6,= -04396.7. b^= 00073.73= 152"-09, 

a, = - -00444,7, a,«' = - -00007,457, 

6, = - -00623,6, 6,*' = - -00010,46 = - 21"-57. 



LECTURE XI. 

THE EQUATION OF THE CENTRE AND THE EVECTION. 

We have seen that the equations of motion 
are satisfied very approximately by the values 



dM dOdl 

d^ dtdt "^ 



V 
I = log- = — Oa COS 2^, 

= n^ + € + 6a sin 2-^, 

where ylr = nt-\- €- (n't + e'), 

and Oa, 62 are small quantities depending upon the ratio nf/n, 
and a is a quantity depending upon n in such a way that 

^, = n« + ^n^'-4 / ''\^ '\'2n'(2n^n')a,\ 
a' 2 32 (w — nf ^ 

while n, € are arbitrary, though subject to the assumption that 
the ratio n/n is small. 

This solution, then, expresses a possible case of motion; 
nevertheless it is no more than a particular case because it 
involves only two arbitrary constants, whereas the complete 
and general solution must contain four, in order that it may 
be able to satisfy any given initial conditions, that is, in order 
that the initial coordinates and their initial velocities may have 
any given values. 
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When there is no disturbance the four arbitrary constants 
are n and e, — which denote quantities similar to those expressed 
by the same symbols above, — and the two elliptic elements e 
and -Bj, of which e denotes the eccentricity of the orbit and 
tsr the longitude of the apse. 

We will now shew how to complete the solution by intro- 
ducing into log(a/r) and additional terms depending on 
quantities similar to e, ^, of which the former is constant and 
the latter varies slowly and uniformly with t ; and for the sake 
of simplicity we will suppose at first that e is so small that its 
square and higher powers may be neglected though it is other- 
wise arbitrary in magnitude. 

Let us assume then 

log - = Og cos 2-^ + e cos (nt — -cj), 

= nt-h€-\-b2 sin 2-^ + 2e (1 + 6©) sin (nt — tsr), 

in which the elliptic terms are of the same form as in the 
undisturbed orbit, and «r is supposed to be slowly variable, 
so that 

where p is supposed to be a small quantity of the order of the 
disturbing force. 

We will now substitute these assumed values in the dilFer- 
ential equations. We have 

-17 = 2 (n — n') tta sin 2-^ + {n'-p)e sin (nt — -bt), 

dH 

-1^ = 4 (ti — nyoq cos 2i|r + (w — pYe cos (nt — tsr), 

J/} 

-17 = n + 2 (n — n') 62 cos 2-^ + 2 (n — /)) (1 + 60) « cos (nt — tsr). 

-1—= — 4 (n — n') 62 sin 2y^ — 2 (n — pY (1 + bo) e sin (nt — nr). 

Hence 
4 (n — ny tta cos 2-^ + (n — pY e cos (nt — -bt) 
+ 2 (n — n') (n—p)ci^ [cos (2-^ — n^ + «r) — cos (2-^ + w^ — 'bj)] 
A. L. 4 
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— {n* + 4rM (w — n') 62 COS 2-^ + 4n (n — p) (1 + 6©) ^ cos (n^ — w) 

+ 4(7l-ll')(w -l')(l+^o)62e[c08(2'<|r-wH t!r)+C0s(2'<|r+n^— or)]} 

+ -3 {1 + 3aa COS 2'<|r} {1 + 3e COS (ne - -or)} 

- n'O 2 + 2 cos 2^ - 3 sin 2>;r [2 (1 4- 60) « sin (n^ - t!r)] ^ = 0, 

and 

- 4 (n - n')« 62 sin 2-^ - 2 (n -p)» (1 + 60) 6 sin (nf - -or) 
+ 4/1 (71 — n!) Og sin 2^ + 2w (ti — p) e sin {nt'-m) 

4-4(n— 7i')(n— p)(l+6o)«a2[sin(2'^— n^+«r)+sin(2'^+?i^-'BT)] 
+ 2 (n — w')(ii — p) efts [— sin (2'<|r — ^^ + -or) + sin (2-^ + n^ — -or)] 

+ n'^ L sin 2>;r + 3 cos 2i^ [2 (1 + 60) e sin (n< - -or)] ■ = 0. 

It will of course be found that with the values of Og, 62 of 
Lecture IV, the terms independent of e vanish identically. 

Equating to zero the coefficients of cos {nt — -or) in the first 
equation and sin {nt — vr) in the second, we get 

(„ _ p)« _ 4n (n -^) (1+ 6.) + ^ = 0, 

-2(w-p)Ml+6o) + 27i(w-p) =0. 

Therefore (n — p) (1 + 60) = n, 

and (n — pf=Am? — ~ 

3 

= u' — ^ n'^y approximately, 

or - = T ^^ = 6o> approximately. 

Now terms have been left outstanding with the arguments 
2-^ — wf + -Bj, 2*1^ + n^ — -BT. These may be removed by assuming 

log - = Oa cos 2-^ + e cos (n^ — 'Bj) + Oai^ cos (2-^ — n< 4- «r) 

+ 022^ cos (2-^ + n^ — -Bj), 

^ = w^ + € + 6a sin 2-^ + 2e (1 + 60) sin {rvb — -bj) 

+ hi^e sin (2i^ — w^ -f- -bj) + 622^ sin (2^ -Vnt — isr). 
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Hence in place of the former equations, we get the 
following 

(«-|,)»-4n(n-p)(l + 6.) + ^ 

Cv 

+ [2 (n - n') (3n - 2«' - p) + 1 g1 a,a„ 

+ \-2(n-n')(n-^'+p)h + ^n''\K 

+ [-2(n-n')(Sn-2n'-p)h + ^n'*\b^=0, 

- 2 (» -pyil + h) + 2n(n-p)- 6n'% (1 + 6,) 

- 2 (n - «') (» - 2n'+p) hiO^ + 2 (m - n') (3« - 2n' - ja) JjOa 

+ [2 (n - n') (n - 2n' +p) a, - 1 n'»1 6a 



+ r_ 2 (n - «') (3» - 2«' - 1>) a, + 1 n'»1 63, = 0. 

Multiply the second by , and add to the first ; this 

•^ w — p 

will eliminate 1 H- 60. 

■^ a' 71— p ^ 

+ 2 (n - n) (n - 2w' +i)) [a2«2i - hK] 
9 « 3 

+ 2 5 [aatta + OaOaa] + g ^'' [&a + 622] 

+ 2 (/l-V)(37l-2w'-p) [02022-63622] 

+ ^ ^"~7; (^ " ^') (^ ~ ^^' + P) [^aOa - a262i] 

- 4 {n - 7i') (3n - 27i' - p) [62O22 - 02622] 



7h 

+ 37l'« [62a-622] = 0. 

71 — p *■ ■* 



4—2 
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Also the equations obtained by equating the coefficients of 
e cos (2-^ — nt + 'Gy) and e sin (2-^ — nt + m) to zero are 

[2 {n - nO (n -p) + 1 ^Ja,-4(n-nO(n-;>)(l +6o)62+3n'Hl+6o) 

+ Kn - 2w' +^)2 + '^1 ^ " 2^ (^ - 2^' + ^) ^21 = 0, 

4 (n - n') (w -p) (1 + 60) ^2 - 2 (n - v!) (n —p) b^ 

-Snf^(l+bo)-(n-2n'+pyb^ + 2n(n-2n' + p)a^ = 0. 

2/1 
Multiply the second by — _a~f > ^^^ ^^^ ^^ *^® first; 

this will eliminate 6a> ^J^d gives 

[2(«-«')(n-p) + |g--Jj:p^(«-r^')(n-i>)(l + 6o)]a. 

Lastly the equations obtained by equating the coefficients of 
e cos (2-^ + nt — tsr) and e sin (2^ + nt — 'ur) 
to zero are 

-2(n-ri0(^-l>)+|5la2-4(n-n')(n-p)(l+6o)62-3w'«(H-6o) 



+ r (3w - 2n' -pf + ^t^ I ^22 - 2n (3n - 2n' ~p) 622 = 0, 



and 



- (3n- 2n' -i)/622 + 2?i(3n - 27i' - p) Oaa = 0. 

2/1 

Multiply the second by — „ ^c-7 and add to the first ; 

^•^ ^ 3n-2n''p ' 

this will eliminate 622* and gives 



XI.] EQUATION OF THE CENTRE AND THE EVECTION. 53 

[_ 2(„ - n')(n -p) + 1 g - 3^^-^(n - nOC" -/»(1 +6.)] a, 
+ [- 4(B - n')(n -p)(l + 6.) - g^^T:^ (» - «')(« - i')] h 

These six equations are to be solved by successive approxi- 
mation ; taking the first rough values of p/n and bo, we find 
from the last two pairs values for On, b^, a^, 622 5 these are 
substituted in the first pair and peld more approximate values 
of p/n and bo, and so on. 

It will be noticed that this complexity is made necessary by 
the &ct that Osq, 621 are found by means of a small divisor 



LECTURE XII. 

THE EVECTION AND THE MOTION OF THE APSE. 

We proceed to the conversion into numbers of the formulae 
of Lecture XI. 



Take n — n'- 


= 1, 


n- 


= 1-08084,9, 


7l' = 


= -08084,9, 




= 1-17150,3, 


log aa = 


= 7-85609, 


log 62 = 


= 800911. 


First Approadmation. 


• 


{n - pf - 


4/1" + 3 ^ = 0, 


4n2 


4-67293,7, 


3': 
0' 


- 3-51450,9, 


{n-pf = 


1-15842,8, 


n—p = 


1-07630,3, 


P = 


-00454,6, 


n — 2n' + p = 


-92369,7, 


Sn-2n'-p = 


3-07630,3, 


l + 6« = 


1 -00422,4 = w/(n-p). 
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Substitute in the equation for a^, 

+6— ?C— (l + 6o) = 0. 
n—2n -\-p^ 

The various terms give 

2 (n - n') (n - p) a, -01545,45 

1 ^. a, 03784,83 

2 a' 

-8 J", (n-n0(n~y)(l + 6o)q« --07264,10 
- 4 (w - nO (n - !>) (1 + 6o) 6» - '04415,05 



n 



^ , — (n-n0(n-l))6a -05144,47 

n - 2n +j9^ ^^ -^"^ 

3n'»(l+6o) -01969,25 



wn'* 



•05373,43 

(ra-2m'+;))> -85321,6 

-4»' + 3^, -116842,8 



- -30621,2 



(% = -17605,6. 
The equation for h^ is 

2?% 7i "■" w 
621 = o / . 021 + 4;; ^r-7- x^(^-p)(l +60)02 

?l — 72 W ^ 
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2n 

Oa -41201,8 



n — 2w' + |) 
n' 



7h ^ ti 

* /„ o„/ ■ „.. (« - P) (1 + 6.) a, 03637,95 



- 2 . " ," ,. (» - p) 6, - 02576,42 



(n-2«'+p)»("--P>^' 



n'» 



- 3 7 sVr-^ (1 + 6.) - -02308,03 

(n — 2w+|)/^ 

6a = -39955,3 
The equation for Oa, is 

r(3n-2n'+p)«-4n«+3^1aa-2(»-n')(«-i))a,+ |^,a, 

I. ^ J ^ 

~^ 3w-2w^-« ^"~"'^^"~^^^^"*"^°^"'~^^"~"'^^"~^^^^"^^''^^'' 
~ * 3w-2n'-j9 ^" ~ "'^ ^^ -P)b,- 3n'' (1 + 6.) 



nn'» 



-«3^3w^a+^.)=o- 

Here - 2 (n - w') (n - p) Oj - •01645,45 

|4«s 03784,83 

2 a' 

~^ 3w-2n'-« ^"~"'^^"~^^^^'^^''^"' --02181,13 

- 4 (n - n') (w -p){l + 6o) 6^ - -04415,05 

- 4 s ^^ (n - n') (n-p)bi - -01544,69 

3n—2n —p 

-3n'»(l + 6„) -01969,25 

- 6 o-^' (1 + W - -01383.78 

- -09254,52 
(3n-2n'-py 946363,4 

-4n»+3^, -1-15842,8 

8-30520,6 
asB= 01114,30. 
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And 



n . M n—n 



* (Sn-i^-py ^"^ - ^> ^^ + ^"^ "» -00327,988 



6a = 01551,37 

Second Approanmation. The complete equation forn — p is 

(„ _p). _ 4n. + 3 ^ + 12 -?^ 6, (1 + 6,) 
'' -^^ a* « — p 

9 u 

+ 2 (» — wO (« - 2n' + p) [OjOji - tgia] + 5 — , [«2«a + «aO»] 

Q 

+ 2^'' [^21 + ^J + 2 (»-riO(3w-2n'-|)) [0^023-62622] 
- 4 {n - n) (3n - 2n' - p) [63022 - 0^622] 

7v ~~' p 

+ 3-^[62i-622] = 0. 

-4w»-f3^. -115842,8 

o' 

12—62(1+60) -00080,8 

n—p ^^ "^ ' 

2 (n - n') (n - 2n' H- p) [03O21 - 62621] - -00520,1 
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I ^3 [0,0^ + 0,0^] -00708,4 

I n'» [6a + 6J -00406,9 



2(n- n') (371 — 2n' — p) [a^o^ — bib^] 



n 



n—p 



(n — n') (n - 2nf + p) [b^a^. - (hbii] - 



- 4 (w — n')(3w — 2nf —p) [b^a^ — 02622] — 

7v """ Xf 



wn'* 



[621 - 622] 



00048,2 
00397,5 

00003,0 
00756,1 



14859,4 



n — p 

(n-pf= 1 

n-._p = 107172,5 

|)= 00912,4 

p : yi = -00844,2. 

Apply these numbers in the equation for 1 + 60 • — 

1+60= 3. To 62(1 +60)-/ ^(^-271' +p)[b^a^-ajb2i] 

n—p {n—py ^ ^ (n—py 

n — n^ 3 n'^ 



n 



n—p 



1-00861,33 



w'a 



-3(^3^2&«(l + ^o) - -00017,51 

^ " ^' (n - 2n' +1)) [ftaOa - Oafca] '00086,55 



(n-p) 
w — n' 



(3n - 2m' - 1)) [ftaOjB - 02622] -00000,64 



(n-py 

3 n'' 
4!(n-py 



[621-622] - -00163,92 



1+60= 100757,1 



Continuing the approximation for a^, 621, 032, 632 the various 
terms found are the following : — 
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•01538,88 
•03784,83 
•07221,53 

•04410,93 

•05097,33 
•01976,81 
•04601,14 

•05365,53 



divisor 
•86169,5 
- 1-15842,8 

- -29673,3 



Oa = 18082,0 



•01538,88 
•03784,83 
-02182,34 

•04410,93 
•01540,41 
•01976,82 
-01390,46 

•09254,01 



divisor 
9-43549,4 
-115842^8 

8-27706,6 



o„ = -01118,03 



•42108,0 
•03598,79 

- 02540,22 

- 02292,94 

b^ = -40873,6 



-00786,805 
-00328,659 
-00231,985 
-00209,403 

6a = -01656,85 



Third Apprommation. We find 



1-15842,8 

- -00081,40 
•00535,02 

- -00723,54 

- -00416,02 
-00048,36 
-00407,40 
-00002,97 

- -00777,57 

(n-py= 1-14838,0 



n-p = 107162.6 

p= -00922,4 

p:n= 00863,6 



-01538,73 
-03784,83 
•07220,54 
•04410,81 
•05096,30 
•01975,96 
-04600,96 

-06365,42 oa = '18092,9. 



divisor 

-86188,04 
- 1-15842,8 
- -29654,8 



-42128,90 
-03597,92 

- -02539,43 

- -02292,60 

6n = -40894,8 



1-00860,75 

- -00017,67 
•00087,94 
•00000,64 

- •OOl 67,84 

1+6,= 100763,9 



- 01538,73 divisor 
•03784,83 9-43488,0 

- -02182,35 -1-15842,8 

- -04410.81 8-27645,2 
-01540,32 

- 01976,96 
-01390,60 

- •09253,93 Oa=^01118,10. 

-00786,881 
-00328,671 
-00231,978 
- 00209,430 

6a =01656,96 
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Fourth Approadmation. 



15842,8 

00081,41 

00534,22 

00727,13 

00416,24 

00048,35 

00407,66 

00002,97 

00778,05 



n — p 

P 

p : n 



1-07160,3 
•00924,6 
•00865,4 



1-00862,79 

- -00017,57 
•00088,04 
•00000,64 

- 00167,9 4 

l+bo= r00766,0 



(n-py= 114833,2 



The values already found for the remaining quantities are 
sufficiently exact. 

These numbers give, taking after Hansen, 

e (1 + 6,) =05491 





e = 


•05449 




2(1 + 


6o)e = 


•10982,0 = 


22651"9 




aae = 


•00986,03 






bae = 


•02228,44 = 


= 4596"-6 




a^e — 


•00060,93 






bae = 


•00084,85 = 


= 175"1, 



and taking the Moon's mean annual motion 17325593", the 
annual motion of the apse is 

148202'' = 4m0' 2". 



LECTURE XIII. 

THE MOTION OF THE APSE, AND THE CHANGE OF 

THE ECCENTRICITY. 

We have seen that when the eccentricity of the Moon's 
orbit is not considered we may write 

l=:l[l+a2Cos2(^-^)], 

3 , 2-m 1 ,3 m* 
where a«=^m'.q .,, ; fh = T-i • 

Let us introduce the two new arbitraries e, «r by writing 
H = hna^[l+hicos(2-2m)0l 

- = Ti-[1 + a2Cos(2 - 2m) ^ + ecos(^-«r)], 

where A is a third arbitrary, which may be chosen to suit our 
convenience ; it must be unity when e = 0. 

Then 
dH dH dt 3 - _ - . Q ,^ /%,\dt 

= - 1 m'TO' ^ [sin (2 - 2«i) ^ + 4«ie cos (2 - 2m) e sin (^ - «)] 

= - 1 m'n' ^'"jl - (4a, + As) cos (2 - 2m) ^ - 4e cos (^ - w)] 
X [sin (2 - 2m) ^ + 4me cos (2 - 2m)^sin(^ - w)]. 
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and also 

^^ = ^ na^ [l-hth co8(2 - 2m)«] -wa%(2 - 2m) h, sin(2 - 2m)0. 

Now we may put A = 1 4- 17, where 17 vanishes with e. Neglecting 
powers of e above the first 

^ = ^=3m^(l+m)esin(r^2^^ + ti7) 
au au 

+ 3m^{l - m) e sin (3 - 2m^ - -or) - 9m^ sin (2 - 2m) ^. 
Neglect at first the last term : 

iy = -3m % J^ ecos(l-2m^+«r)-3m% ~^ eco8(3-2m^-Br). 
1 — 2m o — 2m 

Substitute this in the last term, and we get 

^ = 3m«(l +m)e sin(l-2m^+ m)+ Sm\ 1 -m)esin(3-2m^-'sr) 

CiU 

27 ^ 2 + 4m-4m» . ,^ . 

+ T ^^ (T::2m)T8:r2^) ^ ^"^^^ " "^- 

Now consider the other equation 



— 2m« cos (1 — 2m^ + «?") 4- 2me cos (3 — 2m^ — w)]. 

Differentiate the assumed expression for -, and let h be 

chosen so that the first differential coefficient shall have the 
same form as if h, e, w were constant. 

Thus 
1 dr 1 /o rt \ • /c% -^ \ ^^0 1 ' //I \d0 

where 

- 7- ^ [1 + ^2 cos(2 - 2m)^] + ^ co8(^-«r) + e ^ sin(^-«7)= 0, 
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or 



-^cos(^ — 'cj)4-e ,^ sm(^ — «r) = 6w'(l +m)«sin(l — 2m^ + w) 



and 



+ 6m' ( 1 — m) « sin (3 - 2m0 — «•) 



TT = 7^(2 - 2m) aasm(2 - 2m)^ + ,^ esin (^- «r) 

= ^(2-2m)a,8in(2-2m)^ 

+ -r- [1 + ^ COS (2 - 2m)^ e sin (^ - «), 
^ = -^(2-2m)«a,co8(2-2m)<?^ 

- y (2 - 2m) A, sin (2 - 2ot) ^ e sin (tf - tir) ^ 

+ y[l+A,co8(2-2»»)^ecos(^-«r)^ 

+ ^ [l+A,cos(2-2m)^ |^g8in(d-w)-e^^cos(^-w)l ^. 
Multiply by r^/n^a^ ; then since 



r2~=i7 = Ana«[l+A,co8(2-2m)^], 



we have 



^,^ = (2-2m)'a,cos(2-2m)d 

- (2 - 2»i) ^ sin (2 - 2m) ^ e sin (^ - «■) 

+ e cos (^ — «r) + 2^ cos (2 — 2ot) 9 e cos (0 — w) 

-~(2 - 2m) Oa8in(2- 2m) ^ 

+ [1 +A,co8(2-2m)^'' r^sin(^-Br)-e^cos(^-i«r) I , 
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and this is equal to 



— 2me cos (1 — 2m6 + -cr) + 2me cos (3 — 2m0 - -cr)] 
= l + Y + 2AsCos(2-2m)^|[l+aaCos(2-2m)^+eco8(^+'5T)] 

" "Tl + o ^"^' [1 " 3^2 ^^^ (2 -270)0- Ze cos {6 - w)] 
+ 2 m%M COS (2 - 2m) ^ 4- 6ase COS (^ - "bt) 

- ^2 4- 2m) e cos(l-2m^ + tu) - L - 2m j e cos(3-2m^- tsr) . 

The terms in these two expressions which are independent of e 
give no new information ; equating the others : — 

-(2 - 2m) A^ sin (2 - 2m) ^ sin (^ - m) 
« J^(2 - 2m)a2sin (2 - 2m) ^ 

+ [l4A^cos(2-2m)(9pr^sin(^-tsr)-e^Jco8(^-'cr)1 

3 3 r 

= 3m^ — ^ rri^e cos (^ — id-) + ^ m^ 6a^ cos {0 - id-) 

- f 2 + 2m j e cos (r-"2m^ + «^) - (o " ^m J e cos (3-2m^-tsr) 
+ 9m2i7 cos (2 - 2m) 0. 
Reducing this expression 
- -^ sin (^ - -cr) + e -^ cos (^ - «r) = f 2 m* - jT m* j e cos (^ - «r) 

+ f 5 m^ + 3m3 + ^ mM ecos(l - 2m^ + lar) 
+ ^3m2-3m« + -g-m*)ecos(3-2m^-«r) 
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and from before 
,^ cos (^ — tsr) + e-^ sin (^ — -cr) = lSm*e sin (0 — ^) 

+ (6m'' + 6m0 e sin (F- 2m0 + isr) 



+ (6m2 - 6m») e sin (3 - 2m0 + tsr). 



Hence, 



de / 3 . 147 



= (-|wi« +^m*)esin2(^-i;j) 



27 



+ ( -^m»-3m«- 3 m*^ e sin (2 - 2m) ^ 

/ 15 . 9 3 63 A . ,^ ^ ^ . 
+ ( — -J- m* — X m' — Y^ mM e sm (2m^ — 2«r) 



^ m' - 2 m' - jg mM c sin (4 - 2m0 - 2i!r), 



;e> 



141 . / 3 . 147 



<;(? 4 



= i"^* + ^"^*+( i***' -\|-m')co8 2(d-^) 



+ ( ^m> + 1 m« + ^ m<) cos (2m^ - 2w) 



(Q Q "It \ 

-g^'+o^*"^ 16 ^*)c<^8(4-2m^-2t3-). 

We notice that among these terms, one is of long period, 
approximately semiannual, and will become of greater relative 
importance than the others on integration. 

To eflfect this integration, assume 
n = tu + asin2(^-t;T)-f)8sin(2-2m)^ + 7sin{(4^2m)^-2iar}, 

so that the mean motion of 11 is the same as that of ^, and 
substitute in the equation. 

A. L. 5 
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Then 

X cos 2 (^ — bt) 

r Q 09 3 1 

+ 7 m»- 6w« + -g-m* f(2-2m)/8-6m«a-TmV 

X cos (2 — 2m) 5 
+ r ^m«+|m» + Y|m*-.|m«a1cos(2w5-2n) 
r 3 3 15 9 3 ~l 

X cos {(4-2m)5-2OT}, 
so that if we take 

3 , ^219 , 

^ 9 ,_^15 , 99 , 

3 2 3 8 '51 , 

7= 8^ "16^ 128^' 



we have 
dU 3 , 309 



*+r^m« + |m» + ^m*lcos(2m^~2n). 



d^ 4 ■ 32 
If we write m0 — IT = '^, 

this becomes -j^ = a — 6cos2'^, 

where 

3 3 309 , ,15 .9 3^45 , 

4 oZ 4 Z o 

and the solution is 

tan~^ f . / — , tan •>/r j = 0\/a^ — 6^ + constant. 



1 
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Hence if we denote by -^ the mean rate of change of cr, we have 

da 

3 , 225 , 4071 , 

dwo 3 ,^225 ,4071 , 

T0= 4'"+-32'"+128"'- 

We observe that a + b and a — b are the rates of separation of 
the Sun from the apse when the Sun and the apse are at 
quadratures and syzygies with one another, respectively, — that 
is if we take 11 for the longitude of the apse, or, what is the 
same thing, if we ignore small terms of short period. Hence 
the mean rate of separation of the Sun from the apse is a mean 
proportional between its rates when at quadratures and syzygies 
respectively with the apse*. 

[* This is the analogue for the case of the apse of Machin and Pemberton's 
theorem on the motion of the node, inserted in the third edition of the 
Prineipia as a scholimn to prop, xxziii., lib. in. See some notes by Adams 
in Brewster's Life of Newton^ Appendix zzz.] 



5—2 



LECTURE XIV. 

THE LATITUDE AND THE MOTION OF THE NODE. 

Let us first treat this problem on the supposition that the 
latitude is so small that its square may be neglected. The 
equation of motion, taken from Lecture II, may be written 

where ^^ = r sin (latitude) and the cube of « is omitted; or 
neglecting the parallactic terms 

The value of iijr^ may be considered known by the operations 
which have determined the motion in an orbit coinciding with 
the ecliptic ; that is to say, 

^ = ^ri + |aa2 + 3aaCos2Vr + ^|a2«4 3a4)cos4^|, 

where a has the definition of Lectures IV, V ; or numerically, 
taking 

n — w' = 1, 



And 



^ = 117 150,3 + 02523,0 cos It + 00025,15 cos U. 



^ = n'« = -00653,6. 



Hence 



^ = - ^ [117803.9 + 02523,0 cos 2^ + '00025,15 cos 4^]. 
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Let US now consider the equation 

where P = jo + 2gi cos 2^ -f 2ga cos 4^, 

in which q^, q^ are supposed small. 

Suppose a term in ^r to be ccoB(kt + ^); when this is 
substituted in Pz there will arise terms 



c cos (A? - 2^ + /8) c cos (k + 2t-\- /8) 



c cos (k-U + fi) c cos (A: -f 4^ + /8). 

Let us therefore assume 



z=:c[cos(kt + fi) + Ci cos (A; + 2« + /8) + Ca cos(A? + 4^ + ^) + ... 



+ c_i cos (A; - 2« + )8) + C-3 cos (A; - 4^ + )8) + . . . ] 
c is arbitrary; we have to determine A?, Ci, c_i, &c. 

Substitute and equate coefficients: 

. . . + [-(A:- 4)H?o] 0-2+910-1+92+ • • • =0 

+g'iC_2+[-(A?-2)"+9o]c_i+9i+92Ci + ... =0 

+92C-2 +9iC-i +[-A:'+9o]+9iCi+g'2C2+... =0 

+920-1 + 9i +[-(A?+2)2+9o]ci+9iC2+... =0 

+ 92 +9iCi +[-(A;+4)«+9o]c2+... = 0. 

If 9i, 92,... are neglected, we have simply 

-A;^ + 9o = 0; 
this is a first approximation to the value of k. 

Taking 91 into account and neglecting 92 

c g^ 

'" 9o-(A:-2)»' 

c =-_— ii— 



70 LECTURES ON THE LUNAR THEORY. [lECT. 

In the actual case considered we notice that ^o does not 
diflfer widely from unity. Hence k is nearly equal to unity 
also, and the denominator in Ci is small, and makes d much 
more important than Ci. 

If we substitute these values in the third equation above, 
we have 

^ - ?» + ?'i,T^(kx^ %Tr(^^} = <>' 

whence 

(A;»-?o)'-8(i^-go)'-{16(?o-l) + 29x'}(A»-?o)-8g,» = 0, 
which may be put under the form 

(i^-go? + 2(?o-l)(A:»-9o) = -?i' + j?iHA;^-?o) + |(fc»-(?o)', 
whence 

With this equation we can approximate very rapidly to the 
value of h Taking as a first approximation 

substitute this value of k in the small terms and we get as a 
second approximation 

A: = 1-08516,9. 

Whence the ratio of the retrograde motion of the node to the 
Moon's mean motion is 

A;/n - 1 = 5f - 1 = -00399,7, 

where g is written for kjn. This value is very correct. Taking 
the Moon's mean annual motion as 17325593'^ the resulting 
annual retrograde motion of the node is 

69252" = 19° 14' 12". 
Next find the values of the coefficients c«i, Ci, c_2, Cj. We have 

}i= -01261,5, 
go -(A? -2)2= -34112,3, 
Jo - (A? + 2)« = - 8-34022,8 ; 
whence as a first approximation 

c_i = - -03698,19, Ci = 00151,26. 
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Hence 

?iC-i + 9a = - 00034,02, giCi + ?a = -00014,49, 

and 

go - (A? - 4)» = - 7-31821, Jo - (* + 4)« = - 24*6809, 

whence 

C-, = - 00004,650, 

Ca = 00000,587. 
A second approximation to c_i, Ci gives 

[- (A:- 2y + qo] C-i = - (91 + jaCi + jiC-^), 
[-(A: + 2)2 + 9o]ci = - (91 -f ?3C_i + 9iCa) ; 
with the above values 

9i + 92C1 + 9iC_3 = 01261,47, q^ + 9,0-1 + qi(h = 01261,03, 
so that c_i = - -03698,00, 

Ci = -00151,20. 



LECTURE XV. 

MOTION IN AN ORBIT OF ANY INCLINATION. 

Let us consider the change in the plane of the orbit 
produced in an indefinitely small time dt by the action of a 
given disturbing force. Let Z be the resolved part of the 
disturbing force at any time in a direction perpendicular to 
the plane in which the body is moving at the instant. Imagine 
the force Z to act by impulses at the small intervals of time dt, 
then Zdt will be the indefinitely small velocity generated by 
the force Z in the time dt, in the direction perpendicular to the 
plane of orbit at the instant. 

Let FP be the radius vector and P the position of the 
body at the instant. Also let PT represent the velocity at 




VI 



the instant in magnitude and direction; then if Tt be taken 
perpendicular to the plane FPT and equal to Zdt, the velocity 
and its direction after the impulse will be represented by Pt, 
and the new plane of the orbit by FPt Draw Tm perpen- 
dicular to FP and join tm\ then tmT is the angle through 
which the plane of the orbit has been turned about the radius 
vector FP in the indefinitely short time dt 
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Now 



tmT=^ 



tT Zdt 



Tm 



where v is the resolved part of the velocity at P perpendicular 
to the radius vector. But 

H = vr; 

hence the angle through which the orbit is turned in an 
indefinitely short time dt is 

rZ 



H 



dt 



To find the corresponding changes in the elements that deter- 
mine the plane of the orbit, namely, the inclination of the orbit 
to a fixed plane, and the longitude of the node on that plane. 
Let NPQ be the great circle which represents the plane of the 




orbit at the time t, NR the plane of reference, usually the 
plane of the ecliptic, P the position of the body at the same 
time, and let i = PNRy the inclination, and let N be the 
longitude of the node. 

Let nPq be the position of the orbit at time t + dt. 

Take NQ = 90° ; draw nm perpendicular to NPQ and qQR 
perpendicular to NR\ then Qi2 = i; and by what we have 
just proved 

NPn^^dt 

JlL 

Zr . Zr 

Therefore tvm ^-jjdt, sin ^, qQ = -fTdt, cos 0, 



where 
But 



0=^nP. 
nm = Nn sin i = sin i dN\ qQ = di. 
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Therefore 

di __ Zr C08 6 

It W~' 

dN Zr sine 
dt H sin t ' 

which give the changes of the elements required. 

Now let NM8 be a spherical triangle, the centre of the 
sphere being 0, the centre of gravity of the Earth and Moon ; 




and let GS, OM, ON point respectively to the Sun, the Moon, 
and the node of the Moon's orbit upon the ecliptic, so that NM 
is the plane of the Moon's orbit and N8 the ecliptic. Let 
M8==o), ira= ^, NM= 0y of which the first is identical with 
the quantity denoted by the same symbol in Lecture II, but 
the second and third are not so. 

Then, following Lecture II, the forces on the Moon are 



mr 



7T- 3 cos 0), in SO, 

if we ignore the parallactic terms. 
This latter may be resolved into 

— —^j 3 cos 0) X (cos COS ^' + sin 5 sin ^ cos i) in MO, 

perpendicular to 
^ 3 cos o) X (sin cos 0" — cos ^ sin ^ cos i) MO in the plane 

of the orbit, 

, perpendicular to 

— TT 3 cos a> X sin 0' sin i the plane of the 

orbit. 



mr 
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MOTION IN AN ORBIT OF ANY INCLINATION. 
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Now 

cos 0) = COS COS ^ + sin ^ sin ff cos i 

= C08(5-^)C08»| + C0S(^ + ^8in2|, 

sin 6 cos ^ - cos sin ^ cos i 

= sin (^ - ^) cos' I + sin (^ + ^) sin* I . 

Hence we have the following expressions for the three forces : 

X 2 cos'l sin"! + {1 + cos 2 (^ - ^} sin* I , 



T= 



|^[ 8in2(tf-^)cos«i 



+sin2^.2cos'|sin'|+sin2(^+^)sin*| , 



Z = 



'im'r . . 



-sin(5-2^)cos«^ 



Now we have seen 



+ sin cos i + sin {0 + 2ff) sin= |1 . 



Zr cos 



dN Zram0 



dt H ' dt JBTsini' 

also, the rate of advance of the node along the orbit is 

^sin 
H tan i * 

Thus the equations of motion become 

H^d0 Zr sin 
r» dt if tan i' 



together with 



dt^ r» ' 

dH 



dt 



= -rT. 



LECTURE XVI. 

MOTION IN AN ORBIT OF ANY INCLINATION {continued). 

To satisfy the equation at the end of Lecture XV, assume 

r^a [l+-4iCos 2(^-^)+^2Cos 2^ + ^3 cos 2^+^, cos 2(^+^)], 

neglecting the square of the disturbing force and the eccen- 
tricity ; thus in the small terms we write 



d0 , , dff , 

r^a, -ji=n, r =a, -rr^w- 
at at 



Hence 



- ^ = ri^a [(2 - 2m)Mi cos 2 (^ - ^') + 4^2 cos 25 

+ 4m2^3 cos 20' + (2 + 2myA^ cos 2(6 + 0')] ; 
substitute in the equation 

therefore 

H^=^ fjLa[l -^ A,co82(0 " 0") + A^coa20 + A^co8 20' 

+ A^ cos 2{0 + ^)] 

- 1 mWa" \{1 -f cos 2(0" 0')} cos* | 

+ {cos25+cos2^}2cos2|sin«|+{l-|-cos2(5+^){sin*| 

- nW [(2 - 2m)2 ^1 cos 2 (5 - ^) + 4^2 cos 20 

+ 4^m^As cos 20' + (2 + 2niy A^ cos 2 (5 + 5')]- 
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Again, we have the equation 

dt 
which may be written 

fl- ^ = - 1 n»ni»a« fain 2 (^ - ^ cos* I 

+ sin2^.2sin»|cos«| + sin2(^ + ^)8in*| . 

Substitute for im its approximate value mV in the small terms ; 
and we find from these two equations 

Q ft ^ ?* 

— (1— m)-4i + 4(l — mf -4 1 + ^ m" ( 1 — m) cos* ^ = — ^ m' cos* ^ , 

3 i % 

— ^2 + 4-^2 + o m' . 2 cos^ ^ sin'^ 

3 i . i 

= — ^ m^ . 2 cos* ^ sin^s , 
2 2 2 

3 i i 

- m-43 + 4?m'-48 -f ^ m* . 2 cos^ ^ sin'^ = 0, 

-(l+m)^4 + 4(l+m)»il4 + |mHl+m)sin*^ =-|m«8in*2. 

Therefore 

. _ 3 , ,i 2-m 

^^--^m cos 2^^^^^^^_^^^^-^^_^^^, 

il, = -2m'co8«|8in'|, 

^= am'cos'isin'l^-j-g^lj^-g^^. 



and 
H 



. __3 2 . ,i 2 + m 

*" 2^ ®^^ 2(l+m)(l+2m)(3 + 2m)' 
id 

^« = w'a* Tl + m» - |m» f cos*| + siQ*|) 

+ H ^i COS* H COS 2 (6/ — ^) + 3m cos* ^ sin*^ cos 

21 — m 2 '^ 22 



^li^''-"i~'^<''^*^]- 
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where a is defined by 



fjL = n'a*. 



If we preferred to define a so that the constant term in H^ 
were equal to n^a\ we should have 

3 f i i\ 

nW =^ fjui + rnWa^ — ^ mWa* [cos* ^ + sin* ^ J , 



or 



fi = nW 1 + 5 m" — 3m" sin^ ^ cos' ^ . 



Let us next find the latitude and the motion of the node. 
Suppose that i = t© + At, 

in which At, AjV are small, to is a constant, and N^ varies 
slowly in proportion to the time, so that we may assume 

Ai\r = i\rjsin2(tf-^) + i\r,sin2^ + i\r,sin2^ + JV4sin2(d + ^), 

At =/iC08 2(tf-^) + /,cos2tf + /,cos2^ + /4Cos2(^ + ^'). 
Then remembering that 

cW dN 

an expression that must be aaed in the terms of chief import- 
ance, we have 

^-^=-2(l-m)n/,sin2(tf-^-2n/,8in2^ 

-2(f»-J^W,sin2^-2(l+TO)ii/4sin2(d + ^, 

+ 2 (m-^*^) nA%oo8 2d'+2(l+»»)njy4Cos2(^+^) 

in which tho last tenu will be found to be required to get the 
(u>nHtunt q wrnx^tly to the order m\ 
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These must be equated to 





Zr cos 


Zr sine 




H ' 


Hsini 


respectively. 






Hence 







— g — 2miV,* — mis —r* = — 7 in' cos i ; 

therefore as a first approximatioD 

q = -m^ coat', 
hence 

— 2 (1 - m) /i= Tm^sintcos's » A= — 5 ^ smt cos'^ , 

^ ^ 4 2 81— m 2 

3 3 . . . 

— 2/2 =— -rm'sinicosi, 72= om'sinicosi, 

4 o 



— 2 ( 771+ Tm'cosij /s = — jm'sin i, 



J. ^ 3 m sin I 



8^3 
1 + - m cos I 
4 



(1 + m) 74 = — -7 m'smi sm'^, 74= ^ i si^^* sin'^^ , 

4 2 81+m 2 



and 



2iV2 = ^m'cosi, i\ra= ^wi'cosi, 

0/ . 3 « A ,r 3 . . ,7- 3 mcosi 

2(m+7m'cosiJiV8= ^m^cosi, Ns= ^ ^ , 

1 + T W cos t 

4 

Substitute above for the quantities 7,, -^g and we get the second 
approximation to q. 



5^ = J m' cos i — ^ m' cos' i + ^ m' sin' i 



78 LECTURES ON THE LUNAR THEORY. [LECT. 

where a is defined by 



fi = n^a\ 



If we preferred to define a so that the constant term in H^ 
were equal to n^a\ we should have 



^Q/.4-. 



or 



3 f i i\ 

a* = /^a + rnWa^ — ^ TnWa^ [cos* ^ + sin* ^ J , 

fi = n W 1 + 5 m' — 3m" sin^ - cos" ^ . 



Let us next find the latitude and the motion of the node. 
Suppose that i = t^ + ^i, 

in which Ai, dkN are small, t, is a constant, and Nq varies 
slowly in proportion to the time, so that we may assume 

AiV^ = i\r,sin2(5-^) + J^asin2^ + iVr,sin2^ + JV4sin2(^ + ^), 

Ai =/iCos2(^-^) + /2COs2^ + /sCos2^' + /4Cos2(5 + ^0. 
Then remembering that 

dff dN 

an expression that must be used in the terms of chief import- 
ance, we have 

^ = ^=-2(l-m)n/isin2((?-^0-27i/,sin25 

-2('m-i^)n/38in2^'-2(l+m)n/4sin2(^ + ^), 

^= ^J^ + ^^= - ^+2(l-m)wi\riCos2(^-^)+2ni\r3cos2^ 

+ 2 ^m - - -T^ nN^ cos 2^+ 2(1 + m)nN^ cos 2(5+(?') 

in which the last term will be found to be required to get the 
constant q correctly to the order w?. 
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These must be equated to 
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Zr cos 


ZramO 




~ H ' 


Hsini 


respectively. 






Hence 







-}-2wi\r,2-m/,-^* = -^m"co8i; 
therefore as a first approximation 

gr = -m" cosi; 

hence 

3 . . i r 3 m* . . ,i 
-2(l-m)/i= ^m^sinicos^g, ^i=-81^^®^^* 2' 

3 3 

- 2/2 =-2^^sinicosi, /2= gm^sinicosi, 



— 2 [m+ -m^eos i) I^ = — ^-m^sin i, 



• • 



3 m sin t 



8,3 

1 + Y m cos t 

4j 



-2(l+m)/4 = -|m»smism>|, /<= - j-j-^ sini sin*^ , 
and 

3 3 

2N.^ = 7 m^ cost, i\ra= Qm^cosi, 

/ 3 A TIT 3 . -^ 3 m cos i 

2im-\--^m'^cosijNz- jm^cost, iYs= ^ g , 

^ ^ 1 + T W COS i 

4 

2(l+m)i\r4= ;^m^sin«|, i\r,= -J^sm'^. 

Substitute above for the quantities /s, Ns and we get the second 
approximation to q, 

gr = - m^ cos i — ^rs m^ cos* i + ^r m' sin* i. 
^4 32 64 



80 
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It will be observed that /j, Nl are of lower order than the other 
coeflScients, so that in order to obtain them correctly to the 
same order as the others we were obliged to retain small terms 

in -jT arising from the variability of N. 

If we take the variable plane defined by the longitude of the 
node No and the inclination {q as the plane to which the position 
of the Moon is referred, we have the latitude of the Moon above 
this plane 

= Ai sin — A^ sin i cos 



3 



9 



= Q m sm t cos- ^ 



m 



1 + T w cos i 

4 



1 — m 



sin (tf- 2^0 



— Q m^ sin i cos t sin 

o 



3 . . . « * 
+ Q m sm I sm' - 



m 



- 3 . 1+m 

l-^-rm cos I 
4 



sin {0 + 2ffy 



LECTURE XVIL 
ON hill's method of treating the lunar theory. 

Let us suppose the Moon to move in the plane of the 
ecliptic, and let us refer its motion to rectangular axes in 
rotation, the rotation being such that the axis of x passes 
always through the mean position of the Sun; that is, the 
axes rotate with angular velocity n\ and if we suppose the Sun 
describes a circular orbit about the origin, its coordinates are 

x' = a', y' = 0. 

Let Xy y be the coordinates of the Moon. 

Then the disturbing forces of the Sun upon the Moon 
relative to the Earth are 

m' x^a' m' m' y 

parallel to the axes of x and y respectively, where 

p'^(x^ay + f, 

and the forces of the Earth on the Moon relative to the 
Earth are 





a 


fl X 


7^ r' 


where 




r* = 


a^ + y\ 


Now these forces 


may 


be written 






dfl 


da 






dx ' 


dy' 



where H = - H 



/A m mx 



a'^ 



A. L. 
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p 

Hence 






We have tacitly assumed the origin to be at the centre of the 
Earth ; if we prefer to place it at the centre of gravity of the 
Earth and Moon, the necessary change is effected by multiplying 
the last terms, which correspond to the Parallactic Inequalities, 
by(^-il/)/(i? + if). 

Equating these forces to the accelerations of the Moon 
parallel to the coordinate axes, we have the equations of 
motion in the form 

or, as they may be written, 

dt" dt dx' 

d^y , ft / ^ _ dR 
d^^ di^d^' 

where iJ = fl + ^ ^'"^ (^ + V"") 

Now suppose we have found values of os and y which satisfy 
this pair of equations and which involve two arbitrary constants. 
This may be accomplished by taking assumed developments 

X = ^Oi cos i (t + y), 

y = Xbi sin i {t + 7), 

substituting in the equations, and equating coeflScients of the 
various terms. The solution found will include the Variation 
and the Parallactic Inequalities. Let it be required to ameud 
this solution by the introduction of the remaining two arbitrary 
constants that are required for a complete solution. 



-I- 
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Let the additional tenns that we seek be &r, Sy, which we 
shall suppose so small that their squares and products may be 
neglected, let us consider first the terms which are multiplied 
by the first power of one of the new arbitraries, the original 
particular solution corresponding to the case in which this 
arbitrary is zero. 

Then Bx, 8y are determined by the equations 

where X, Y are supposed known functions of x, y or of t, and 

have been added here to include disturbing causes not allowed 

for in the above form of R, 

dx ui/ 
Multiply the original equations ^Y -Jiy -ri and add : 

d?x dx d^y dy _ dR dx dR dy 
d^diWdi''dxdi'^dydi - 

_dR 

" dt ' 

since x, y are the only functions of t that R involves ; whence 



(sy-d)-^-". 



where C is an arbitrary constant ; this is the integral known as 
Jacobi's Integral. 

Let us write 

then we have V' = 2R + C; 

and from the original equations themselves 
dV d^x . ,d?y . , 

dR . dR . , 

= ^- cos + ^- sm ; 

dx dy ^ 1 

i 
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and 

umXj d'u 
And from these, differentiating and substituting for -^ , -^ , 

we get 

d'V 
dt^ 



-^^{^*^'')-^[ s«-* 












LECTURE XVIII. 

ON hill's method of treating the lunar 

THEORY (continued). 

The equations for Sw, hy are 

dt* dt da^ dxdy ^ 

d^hy „ ,dhc d^R ^ d^R ^ ,^ 
W^^^-dt^l^y^^^-dy^ ^y^^'^ 

the equations for Xy y are 

d^" di^dco' 

dt'^ dt" dy' 

dm (lAJ 
Multiply the former pair by -^, -# respectively, and the 

latter pair by -i— , -^, and add all together; we get 

dxd^ fo dhx dy d^'By d^y dhy 
dt dt^ "^ dt^ dt dt dt^ ^W dt 

^(dJ'Rdx d^ dy\^ fd^dx d^R dy\ ^ 
" Kdx' dt ^ dxdy dt) '^ \dxdy dt "^ dy^ dt) ^ 

dR dhx dR dhy y dx ^dy 
dx dt dy dt dt dt ' 

No ^ ^ + ^L^ dj_d dR 

da^ dt dxdy dt " dt dx' 

d^R dx d^Rdy _d dR 
dxdy dt dy^ dt "" dt dy ' 

- - --- :-' :--- 
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Then our equation may be integrated 

da; dSx dy dSy _dR^ dR^ ,p 
diW^S dt~ dx^^ dy^y^^' 

SO that T is a known function of t, which involves an arbitrary 
constant. 

Now let us assume 

hx = v cos ^ — ii; sin <^, 
Sy = v sin (fy-k-w cos (f>. 

Substitute above for -n 7 -j. , -rr , -jf- ', we find 
yrfdv dd>\ fdR , dR , ,\ 

/ dR . , dR ,\ rn 

+ I — -1- 8in<pH- 7- cos 0) w + i. 

^ , dR , dR . , dV 

But _co8^ + -^sin^ = -^, 

dR . ^ dR , Trfd6 « A 
Therefor F(*-.f).f,+ r(f + 2»')».r, 

whence ^= JV\^ "^ ^') ^^^ "^ I Fi ^^— 

An arbitrary constant is included on the right. This 
equation shews that when w is known, v can be found ; it 
remains to determine w. 

Now by actual differentiation 
. dBx . . . dBy dv dd> 



dt ^ dt dt dt 

d^Bx . __ ^ d% d^w . ^dv d0 /c^<f)^2 . cP^ 



-9m0^ + cos<A^/=-^^,-f.2^^-^(^)+. 



